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PREFACE 


1.  The  Corps  of  Engineers,  through  its  Civil  Works  program, 
has  sponsored,  over  the  past  23  years,  research  into  the  behavior 
and  characteristics  of  tidal  inlets.  The  Corps'  interest  in  tidal 
inlet  research  stems  from  its  responsibilities  for  navigation, 
beach  erosion  prevention  and  control,  and  flood  control.  Tasked 
with  the  creation  and  maintenance  of  navigable  U.  S.  waterways, 
the  Corps  routinely  dredges  millions  of  cubic  yards  of  material 
each  year  from  tidal  inlets  that  connect  the  ocean  with  bays, 
estuaries,  and  lagoons.  Design  and  construction  of  navigation 
improvements  to  existing  tidal  inlets  are  an  important  part  of 
the  work  of  many  Corps  offices.  In  some  cases,  design  and  con- 
struction of  new  inlets  are  required.  Development  of  information 
concerning  the  hydraulic  characteristics  of  inlets  is  important 
not  only  for  navigation  and  inlet  stability,  but  also  because 
inlets,  by  allowing  for  the  ingress  of  storm  surges  and  egress 

of  flood  waters,  play  an  important  role  in  the  flushing  of  bays 
and  lagoons. 

2.  A research  program,  General  Investigation  of  Tidal 

Inlets  (GITI),  was  developed  to  provide  quantitative  data  for 

use  in  design  of  inlets  and  inlet  improvements.  It  is  designed 

to  meet  the  following  objectives: 

To  determine  the  effects  of  wave  action,  tidal 
flow,  and  related  phenomena  on  inlet  stability 
and  on  the  hydraulic,  geometric,  and  sedimentary 
characteristics  of  tidal  inlets;  to  develop  the 
knowledge  necessary  to  design  effective  naviga- 
tion improvements,  new  inlets,  and  sand  transfer 
systems  at  existing  tidal  inlets;  to  evaluate 
the  water  transfer  and  flushing  capability  of 
tidal  inlets;  and  to  define  the  processes  con- 
trolling inlet  stability. 

3.  The  GITI  is  divided  into  three  major  study  areas: 

(a)  inlet  classification,  (b)  inlet  hydraulics,  and  (c)  inlet 
dynamics . 
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a.  The  objectives  of  the  inlet  classification  study  are 
to  classify  inlets  according  to  their  geometry, 
hydraulics,  and  stability,  and  to  determine  the 
relationships  that  exist  among  the  geometric  and 
dynamic  characteristics  and  the  environmental 
factors  that  control  these  characteristics.  The 
classification  study  keeps  the  general  investiga- 
tion closely  related  to  real  inlets  and  produces  an 
important  inlet  data  base  useful  in  documenting  the 
characteristics  of  inlets. 

b.  The  objectives  of  the  inlet  hydraulics  study  are  to 
define  the  tide-generated  flow  regime  and  water- 
level  fluctuations  in  the  vicinity  of  coastal  inlets 
and  to  develop  techniques  for  predicting  these 
phenomena.  The  inlet  hydraulics  study  is  divided 
into  three  areas:  idealized  inlet  model  study, 
evaluation  of  state-of-the-art  physical  and  numer- 
ical models,  and  prototype  inlet  hydraulics. 

(1)  The  idealized  inlet  model.  The  objectives  of 
this  model  study  are  to  determine  the  effect 
of  inlet  configurations  and  structures  on  dis- 
charge, head  loss,  and  velocity  distribution 
for  a number  of  realistic  inlet  shapes  and  tide 
conditions.  An  initial  set  of  tests  in  a 
trapezoidal  inlet  was  conducted  between  1967 
and  1970.  However,  in  order  that  subsequent 
inlet  models  are  more  representative  of  real 
inlets,  a number  of  "idealized"  models  repre- 
senting various  inlet  morphological  classes  are 
being  developed  and  tested.  The  effects  of 
jetties  and  wave  action  on  the  hydraulics  are 
included  in  the  study. 

(2)  Evaluation  of  state-of-the-art  modeling  tech- 
niques. The  objectives  of  this  part  of  the 
inlet  hydraulics  study  are  to  determine  the 
usefulness  and  reliability  of  existing  physical 
and  numerical  modeling  techniques  in  predicting 
the  hydraulic  characteristics  of  inlet/bay 
systems,  and  to  determine  whether  simple  tests, 
performed  rapidly  and  economically,  are  useful 
in  the  evaluation  of  proposed  inlet  improve- 
ments. Masonboro  Inlet,  North  Carolina,  was 
selected  as  the  prototype  inlet  which  would  be 
used  along  with  hydraulic  and  numerical  models 
in  the  evaluation  of  existing  techniques.  In 
September  1969  a complete  set  of  hydraulic  and 
bathymetric  data  was  collected  at  Masonboro 
Inlet.  Construction  of  the  fixed-bed  physical 
model  was  initiated  in  1969,  and  extensive 
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tests  have  been  performed  since  then.  In  addi- 
tion, three  existing  numerical  models  were  ap- 
plied to  predict  the  inlet's  hydraulics.  Ex- 
tensive field  data  were  collected  at  Masonboro 
Inlet  in  August  1974  for  use  in  evaluating  the 
capabilities  of  the  physical  and  numerical 
models . 

(3)  Prototype  inlet  hydraulics.  Field  studies  at  a 
number  of  inlets  are  providing  information  on 
prototype  inlet/bay  tidal  hydraulic  relation- 
ships and  the  effects  of  friction,  waves,  tides, 
and  inlet  morphology  on  these  relationships. 

c.  The  basic  objective  of  the  inlet  dynamics  study  is 
to  investigate  the  interactions  of  tidal  flow,  inlet 
configuration,  and  wave  action  at  tidal  inlets  as  a 
guide  to  improvement  of  inlet  channels  and  nearby 
shore  protection  works.  The  study  is  subdivided 
into  four  specific  areas:  model  materials  evalua- 
tion, movable-bed  modeling  evaluation,  reanalysis 
of  a previous  inlet  model  study,  and  prototype  in- 
let studies. 

(1)  Model  materials  evaluation.  This  evaluation 
was  initiated  in  1969  to  provide  data  on  the 
response  of  movable-bed  model  materials  to 
waves  and  flow  to  allow  selection  of  the  opti- 
mum bed  materials  for  inlet  models. 

(2)  Movable-bed  model  evaluation.  The  objective 
of  this  study  is  to  evaluate  the  state-of-the- 
art  modeling  techniques,  in  this  case  movable- 
bed  inlet  modeling.  Since,  in  many  cases, 
movable-bed  modeling  is  the  only  tool  available 
for  predicting  the  response  of  an  inlet  to  im- 
provements, the  capabilities  and  limitations  of 
these  models  must  be  established. 

(3)  Reanalysis  of  an  earlier  inlet  model  study.  In 
1957,  a report  entitled,  "Preliminary  Report: 
Laboratory  Study  of  the  Effect  of  an  Uncon- 
trolled Inlet  on  the  Adjacent  Beaches,”  was 
published  by  the  Beach  Erosion  Board  (now  CERC) . 
A reanalysis  of  the  original  data  is  being 
performed  to  aid  in  planning  of  additional 

GITI  efforts. 

(4)  Prototype  dynamics.  Field  and  office  studies 
of  a number  of  inlets  are  providing  information 
on  the  effects  of  physical  forces  and  artificial 
improvements  on  inlet  morphology.  Of  particular 
importance  are  studies  to  define  the  mechanisms 
of  natural  sand  bypassing  at  inlets,  the 
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response  of  inlet  navigation  channels  to 
dredging  and  natural  forces,  and  the  effects 
of  inlets  on  adjacent  beaches. 

4.  This  appendix  discusses  the  calibration,  base  tests, 
and  predictive  test  of  a numerical  model  of  Masonboro  Inlet, 

N.  C. , conducted  as  part  of  the  evaluation  of  the  state-of-the- 
art  inlet  modeling  techniques.  It  presents  the  data  necessary  for 
a comparison  of  results  of  the  physical  and  numerical  models  dis- 
cussed in  the  basic  report  and  in  the  following  appendixes: 

a.  Appendix  1.  R.  A.  Sager  and  W.  C.  Seabergh,  "Fixed- 
Bed  Hydraulic  Model  Results." 

b.  Appendix  3.  R.  M.  Chen  and  L.  A.  Hembree,  "Numer- 
ical Simulation  of  Hydrodynamics  (TRACOR)." 

c.  Appendix  4.  C.  J.  Huval  and  G.  L.  Wintergerst, 
"Simplified  Numerical  (Lumped  Parameter)  Simulation." 
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ABSTRACT 


A two-dimensional  finite  difference  tidal  hydrodynamics 
model  has  been  adapted  to  Masonboro  Inlet,  North  Carolina.  The 
model  was  Initially  calibrated  to  simulate  prototype  conditions 
of  September  and  October  1969.  Operational  runs  were  then  under- 
taken to  simulate  tides  and  currents  corresponding  to  pre-project 
undeveloped  Inlet  conditions  of  November  1964  and  modified  Inlet 
and  north  jetty  conditions  of  June  1967.  For  each  case,  the 
model  was  operated  with  mean  and  spring  ocean  tides. 

The  tidal  hydrodynamics  model  applied  at  Masonboro  Inlet 
uses  an  explicit  numerical  solution  of  the  basic  equations  of 
motion  and  the  continuity  equation.  As  structured,  the  model 
includes  bathymetric  data  for  the  Inlet,  bottom  roughness, 
inundated  areas  and  other  features  of  the  system.  Because  of 
the  size  and  detail  required  In  the  simulations,  a dual-model 
approach  Is  employed.  In  this  technique  a coarse  grid  model 
Is  utilized  to  compute  tidal  flows  which  In  turn  are  applied 
as  boundary  conditions  to  a finer  sub-model  of  the  Immediate 
Inlet  area.  Special  requirements  necessary  for  sequential  oper- 
ation of  the  two  models  are  presented. 
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I.  INTRODUCTION 


Background 

This  study  was  initiated  by  the  U.  S.  Army  Coastal  Engineering  Research 
Center  to  help  evaluate  the  degree  to  which  mathematical  models  can  be  used 
to  predict  quantitatively  the  hydrodynamics  of  flow  through  tidal  inlets 
(exclusive  of  sediment  transport).  The  project  is  a basic  part  of  a more 
extensive  investigation  to  determine  the  effectiveness  of  both  mathematical 
and  physical  models  for  predicting  tidal  amplitudes  and  current  velocities 
for  selected  prototype  conditions.  As  a "typical"  problem  area,  Masonboro 
Inlet,  North  Carolina,  for  which  some  actual  tide  and  flow  data  are  available, 
has  been  selected  for  comprehensive  field  and  model  investigations. 

A specific  requirement  of  this  work  was  that  prescribed  tidal 
hydrodynamic  simulations  (including  calibrations  and  applications)  be  under- 
taken using  existing  mathematical  modeling  procedures.  For  this  purpose, 
HYDTID , a two-dimensional  finite  difference  computational  model  was  applied 
to  Masonboro  Inlet.  HYDTID,  with  its  genesis  in  the  hurricane  surge  model 
of  Reid  and  Bodine  [1968],  was  formulated  and  programmed  as  a basic  part  of 
a comprehensive  study  for  the  development  of  estuarine  transport  models, 

Masch,  et  al  [1969]  and  Masch  and  Brandes  [1971],  and  has  been  developed  to 
its  present  form  through  a series  of  application-improvement  efforts.  The 
basic  outputs  of  this  model  are  the  time  histories  of  tidal  amplitudes  and 
vertically  integrated  velocities  (or  flows  per  unit  width)  in  each  of  two 
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areawise  coordinate  directions.  The  model  utilizes  as  inputs  excitation  tides, 
inflows  and  wind  stress,  and  accounts  for  irregular  bathymetry,  bottom  rough- 
ness, islands,  inundated  areas  and  other  physiographic  features.  Utilization 
of  HYDTID  in  a variety  of  prototype  applications  by  Brandes  and  Masch  [1971a, 
1971b,  and  1972]  has  resulted  in  an  expansion  of  the  model's  computational 
and  output  capabilities  and  improved  programming  efficiency. 


Objectives  and  Scope 


The  primary  objective  of  this  project  was  to  apply  the  two-dimensional 
tidal  hydrodynamic  modeling  capabilities  contained  in  HYDTID  to  Masonboro 
Inlet,  North  Carolina.  More  specifically,  the  project  provided  for: 

1.  Adaptation  of  HYDTID  to  the  geometric  and  bathymetric  configuration 
of  Maeoriboro  Inlet  in  October  1969. 

2.  Calibration  of  HYDTID  for  prototype  tides  and  currents  of  12  Septem- 
ber 1969. 

3.  Application  of  the  calibrated  model  to  predict  tides,  currents  and 
tidal  prisms  at  selected  stations  for: 

a.  Pre-project  undeveloped  inlet  conditions  corresponding  to 
November  1964  data,  and 

b.  Modified  inlet  and  north  jetty  conditions  corresponding  to 
June  1967  survey  data. 

In  each  application,  the  model  was  operated  with  prescribed  mean 
and  spring  ocean  tides. 

Additionally  in  this  report  the  modeling  capabilities  in  HYDTID  are 
described  and  details  of  the  basic  equations,  boundary  conditions,  numerical 
solution  scheme  and  programming  techniques  are  presented.  This  is  followed 
by  the  application  of  HYDTID  together  with  a discussion  of  the  requirements 
imposed  by  Masonboro  Inlet.  Special  attention  is  given  to  the  resolution, 
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time  steps  and  requirements  for  computer  time  and  storage. 

The  HYDTID  model  as  applied  to  Masonboro  Inlet  is  programmed  in  FORTRAN 
IV  source  language  at  a level  compatible  with  that  on  the  UNIVAC  1108  computer 
system.  The  operating  characteristics,  the  I/O  data  structures,  support  data 
management  procedures  and  output  report  structures  necessary  to  implement  the 
model  and  evaluate  the  results  are  provided  in  the  Appendix. 

Description  of  Study  Area 

Masonboro  Inlet  and  the  surrounding  area  to  which  the  model  has  been 
applied  is  shown  in  Figure  1.  The  topographic,  physiographic  and  hydro- 
graphic  details  are  contained  in  U.  S.  Geological  Survey  Topographic  Map 
No.  N3400-W7745/15  and  U.  S.  Coast  and  Geodetic  Survey  Hydrographic  Chart 
No.  1235.  Supplemental  information  including  hydrographic  surveys  for  Novem- 
ber 1964,  June  1967,  and  October  1969,  and  pertinent  details  on  the  north 
jetty  and  its  overflow  section  were  provided  by  the  Coastal  Engineering 
Research  Center. 

The  immediate  study  area  shown  in  Figure  1 covers  approximately  nine 
square  miles  including  an  ocean  approach  section.  It  is  characterized  by  a 
combination  of  inter-connected  channels  including  the  primary  inlet,  Masonboro 
Channel,  Banks  Channel,  Shin  Creek  and  the  Atlantic  Intracoastal  Waterway. 

It  is  also  seen  that  considerable  area  inshore  of  the  inlet  is  subject  to 
Inundation  above  mean  low  water  (MLW).  This  causes  flows  in  some  locations 
to  be  generally  confined  to  the  identified  channels  when  tides  are  near  MLW, 
whereas  when  tides  are  near  mean  high  water  (MHW),  large  areas  are  inundated 
and  flows  are  less  well  defined. 
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Figure  I 


LOCATION  MAP  OF  MASONBORO  INLET  AND  STUDY  AREA 
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Although  not  shown  in  Figure  1,  the  study  area  probably  receives  some 
local  inland  drainage  and  return  flows  through  Bradley  and  Hewlet  Creeks. 
However,  no  data  on  the  magnitudes  of  these  flows  were  available,  nor  is 
there  any  information  on  their  distribution  in  the  inland  waterway  and  out 
through  Masonboro  Inlet. 
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II.  FORMULATION  AND  SOLUTION  TECHNIQUE 


Tidal  Hydrodynamic  Equations 

Mathematical  characterization  of  the  hydrodynamics  in  a tidal  system 
requires  the  simultaneous  solution  of  the  dynamic  equations  of  motion  and  the 
unsteady  continuity  equation.  The  theoretical  basis  for  these  equations  has 
been  dealt  with  in  detail  by  Dronkers  [1964]  and  will  not  be  repeated  here. 

If  it  is  assumed  that  vertical  accelerations  are  negligible,  pressures 
are  hydrostatic  over  the  depth  and  the  fluid  density  is  homogeneous,  the 
equations  of  motion  for  a two-dimensional  (areawise)  tidal  flow  system  can 
be  written  as: 


9t  9x  d 9y 


cos  V 


qx  3q 

d~  9^ 


q^  9q^ 

d 9y 


+ fiq* 


gd  § - gdSey  + K Vw  sin  y (2) 


The  equation  of  continuity  for  the  corresponding  flow  system  can  be  expressed 
as: 


In  eqs.  (1),  (2),  and  (3),  q and  q are  the  vertically  integrated  flows  per 

a y 

foot  of  width  at  time  t in  the  x and  y directions,  respectively  (x  and  y 
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taken  in  the  plane  of  the  surface  area);  h is  the  water  surface  elevation  with 
respect  to  mean  sea  level  (MSL)  as  datum;  d is  the  depth  of  water  at  (x,  y,  t) 
and  is  equal  to  (h-z)  where  z is  the  bottom  elevation  with  respect  to  MSL  meas- 
ured negatively  downward;  r is  the  rainfall  intensity;  e is  the  evaporation 
rate;  Vw  is  the  wind  speed  at  a specified  elevation  above  the  water  surface; 

V is  the  angle  between  the  wind  velocity  vector  and  the  x-axis;  and  K is  a non- 
dimensional  average  wind  stress  coefficient  taken  as  1 .8  x 10~6,  a value  consist- 
ent with  those  suggested  by  Roll  [1965]  and  Wilson  [1960]  for  light  winds.  If 
it  is  further  assumed  in  eqs.  (1)  and  (2)  that  the  dissipation  of  momentum  due 
to  bottom  friction  and  turbulence  can  be  expressed  by  the  Manning  Equation,  then 
Sex  = (n2/2-21  dl0/3)  9 * 9X  and  Sgy  = (nz/2 .21  d10/3)  q • qy  where 

q = (q*  + q^)1/2  and  n is  the  Manning  roughness  coefficient  estimated  from 
x y 

formulas,  tables  or  comparisons  of  measured  and  computed  tide  and  velocity 
variations.  The  Coriolis  parameter,  ft,  is  equal  to  2 u>  sin  <t>  where  w is  the 
earth's  angular  velocity  taken  as  0.73  x 10_u  radians/second  and  $ is  the 
latitude  equal  to  34.2°  at  Masonboro  Inlet. 

Solution  Scheme 

The  basic  tidal  hydrodynamic  equations  presented  above  are  nonlinear 

partial  differential  equations  in  which  there  is  also  a dependence  of  d and 

q on  the  values  of  the  three  unknowns,  q , q , and  h.  This  difficulty,  to- 

x y 

gether  with  the  complex  geometry,  various  interior  features  and  variable 
boundary  conditions  found  in  actual  tidal  systems,  makes  purely  analytical 
approaches  impractical.  For  these  reasons,  numerical  schemes  implemented  on 
high  speed,  large  capacity  digital  computers  can  be  used  to  obtain  approxi- 
mate solutions  to  eqs.  (1),  (2),  and  (3). 
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The  primary  advantage  of  the  numerical  approach  is  that  it  permits 
discretization  of  a prototype  system  and  facilitates  prescribing  in  a more 
accurate  manner  the  inputs  necessary  to  describe  the  hydrodynamic  behavior. 

A numerical  model  consists  of  a series  of  finite  elements  arranged  in  time 
and  space  so  that  the  output  from  one  element  becomes  the  input  to  the  next 
and  so  on.  Each  input  to  an  element  is  operated  on  by  a transfer  function, 
and  through  an  advancing  series  of  spatial  and  time  steps,  the  functional 
behavior  of  an  entire  prototype  can  be  approximately  determined.  Thus,  it 
becomes  possible  to  evaluate  the  overall  response  of  a total  system  by  con- 
sidering the  interaction  between  individual  elements  satisfying  common  boundary 
conditions  in  succession. 

It  is  important  to  recognize  that  the  precision  of  the  results  obtained 
by  numerical  procedures  depends  in  part  upon  the  time  interval  and  the  element 
size  selected  for  the  calculations  and  the  rate  of  change  of  the  phenomena 
being  studied.  The  greater  the  number  of  finite  time  intervals  used  over  the 
total  period  of  investigation  and  the  greater  the  number  of  spatial  elements 
used  to  describe  physical  detail,  the  greater  also  will  be  the  accuracy  of 
the  expected  results.  Thus,  solutions  obtained  may  be  refined  within  the 
basic  assumptions  to  meet  the  demands  of  accuracy  by  decreasing  the  time 
interval  and  element  size,  although  such  reductions  can  be  made  only  at  the 
expense  of  increased  operational  costs.  In  addition  to  constraints  imposed 
on  the  solution  method  by  economics  or  by  desired  accuracy,  there  are  also 
mathematical  considerations  involving  stability,  convergence  and  compatibility 
which  limit  the  computational  time  interval  and  element  size.  It  is  only 
when  a satisfactory  compromise  of  all  these  factors  is  achieved  that  a potenti- 
ally useful  simulation  tool  results. 
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The  computational  scheme  used  in  the  model  to  solve  the  two  equations 

of  motion  and  the  unsteady  continuity  equation  involves  a straight  explicit 

formulation.  In  this  method,  eqs.  (1),  (2),  and  (3)  are  written  in  finite 

difference  form  and  the  three  basic  unknowns,  q , q and  h,  are  determined 

* y 

for  each  computational  element  at  an  advanced  time  level  in  terms  of  known 
conditions  at  the  present  time  level.  This  solution  scheme  is  similar  to  that 
used  by  Reid  and  Bodine  [1968]  in  their  hurricane  surge  model  of  Galveston 
Bay. 

To  facilitate  application  of  this  solution  method,  a computational  grid 
composed  of  a mesh  of  square  elements  in  two-dimensional  Cartesian  coordinates 
is  used  to  represent  the  physical  characteristics  of  the  prototype  system. 

The  solution  of  the  unknowns  for  each  individual  grid  element  proceeds  spatially 
in  a row  by  row  manner  over  the  entire  array  for  a prescribed  set  of  boundary 
conditions.  Time  is  then  advanced  and  the  entire  process  repeated  for  a new 
set  of  initial  and  boundary  conditions. 

Variable  definitions  used  in  this  explicit  method  are  illustrated  on  the 
discrete  element  in  Figure  2.  Each  element  of  this  type  is  identified  by  the 
indices  (i,j),  i representing  the  x-direction,  and  j the  y-direction,  both 
coordinates  being  in  the  horizontal  plane.  The  indices  i and  j increase 
with  positive  x and  y,  respectively.  In  the  numerical  analog,  the  discharges 
per  unit  width  in  the  x and  y directions  are  defined  at  the  centers  of  the 
right  and  upper  sides  of  each  cell,  respectively.  The  value  of  h is  taken  as 
the  mean  sea  level  water  surface  elevation  and  is  defined  at  the  center  of  the 
cell.  Also  defined  at  the  center  of  each  cell  are  the  bottom  elevation,  the 
water  depth,  the  Manning  roughness  coefficient,  and  the  rainfall  and  evaporation 
rates. 


9 


qy(i,j+D 


q«(i,j+D 


qy(i+i,J) 


* 


d(i+  l,j) 
h(i+l,j) 
z(i+  l,j) 

•(i+l,j)l 

n(i  + l,j) 


q«(i  + l,j)-H 


CELL(i  + l,j) 


Figure  2 


DISCRETE  ELEMENT  VARIABLE  L0CATI0N6  USED  IN  'HYDTID' 
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The  explicit  solution  method  used  is  a time-centered  difference  scheme 

involving  time  operations  of  the  "leap  frog"  type  for  computation  of  flows 

and  water  levels.  For  illustrative  purposes,  the  following  time  notations 

are  used  to  develop  the  recursion  relations:  t-1  = (k-1/2)  At;  t - K At; 

t+1  = (k+1/2)  At;  and  t+2  = (k+1)  At  where  k is  an  integer.  Time  staggered 

computations  are  made  for  q and  q at  odd  (k+1/2)  time  levels  and  for  h at 

x y 

even  (k,  k+1)  time  levels.  The  wind  stresses  are  applied  at  spatial  locations 

consistent  with  q and  q but  at  even  time  levels.  The  selected  time  differ- 
x y 

encing  scheme  is  such  that  the  differences  (q*+1  - q^"1 ) and  (q^+1  - q^’1^ 

xx  y y 

are  centered  in  time  at  the  level  of  h*  and  the  difference  (ht+^  - h*)  is 

centered  in  time  at  the  level  of  q^+1  or  qj'+"' . 

x y 

Basic  Finite  Difference  Equations 


Using  the  above  notations  and  writing  the  derivatives  as  centered 
differences,  the  terms  in  the  equation  of  motion  for  the  x direction,  eq. 
(1),  can  be  approximated  as  follows: 


3qx  _ qj+1(i.j)  - qj_1(i.j) 

at  = At 


(4) 


qx  3qx  q*+1(i.j)  q^_1(i+l.  j)  - q^“1(i-1,  j) 

(F  3x“  ” cT  2 Ax 

where  dx  = A™,  J) 

qy  9qx  _ qy  qxt'1(i-  j+1)  ' qx"1(1, 

d dy~  " 3^  2 Ay 


(5) 


(6) 
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where  qy  = j-[qj_1(i,j)  + qj_1(i+l,j)  + qj-1(i,j-l) 


+ qy_1(i+l,  j-1)] 


n H 3h 

9 d 3x 


• !1  q 


. g d ^1111  f JJL.  ~ 

y v AX 


gdS 


ex 


9 n2  qx  q/d 
2.21  d4/3 


q nv 


2.21  dx4/3  x 


(7) 


(8) 


[q!+1(i,j)  • q7B]  (9) 


ni  i + ni+l  i 
where  nx  = — ^ ^ -2-  * 


q7d 


(q!_1(i,j)2  + q2)l/2 


In  these  expanded  terms,  the  only  unknown  quantity  is  qx+^ , the  flow 
per  unit  width  in  the  x direction  at  time  level  (k+1/2)  At.  Combining  the 
terms  according  to  eq.  (1)  and  solving  for  qx  \ a direct  explicit  expres- 
sion for  qx  can  be  written  as  follows: 

t+i,  , t-i  9v  q!-1^*  J+1)  - q!_1(i.  j-D 
qj1(i.j)  = [qj'(1.j)  - Atg^  x 


+ At  fi  qy  - At  g _x 


2 Ay 

d h^l,  j)  - Vli.j) 


Ax 


+ At  KM2  cos  0]  /C 


fx 


(10) 
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where,  Cfx 


1 + 


4t 

T%  TS 


+ 


9 nx  At 

2.21  dx4/3 


CqTdT] 


01) 


A similar  expression  for  qy  at  time  level  (k+1/2)  At  can  be  developed 
from  the  equation  of  motion  in  the  y direction,  eq.  (2),  and  written  in 
difference  form  as  follows: 


qj+1  (i  »j ) = [q^"'(i,j)  - At 


.t-l„ 


qv  q!_1(i+l,  j)  - q^d-l,  j) 


2 Ax 


g d - ^d.j) 

a y 

+ At  K W2  SIN  0]  /C, 


- At«qx  -At,-y  ^ 


'fy 


At  qj"  (1,  j+1)  - qj‘  (1.  j-1) 

-1""  cf,  ■ ' * ^ * rrr1 


9 % At 

+ — * — COT 

2.21  dy4/3 


dt(i «j)  + dt(i,j+l) 


qx  = T^-1(1,j)  + q^O.j+l)  + qj_1  (i-1 , j) 


+ qj_1(i-l,  J+l)  ] 


(12) 


(13) 


13 


ni  ,J  * ni , 3+1 


q7d 


(q!  + q!_1(1,j)2)l/2 


Similarly,  the  third  unknown,  the  water  level  h,  can  be  determined  from 
the  continuity  equation,  eq.  (3)  by  approximating  its  derivatives  with  centered 
differences  as  follows: 


3h  _ ht+2(i,j)  - ht(1,.1) 

at  “ At 


3qx  _ qx+1(i,j)  - qx+1(i-i,  j) 

ax  " Ax 


3qy  _ qr1(i,j)  " ^+1{i*j_1) 

ay  Ay 


(14) 


(15) 


(16) 


t+2 

Combining  these  difference  terms  and  solving  for  h , the  following 
expression  for  water  level  results, 


ht+2(i,j) 


t qx+1(i,j)  - qy+1(i_1»  j) 

h il  ,j)  - At  


Ax 


qj+1(i,j)  - q*+1(1,  j-1) 
- At  . J- 


Ay 


+ (r  - e)  At  (17) 
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It  should  be  noted  that  in  developing  eq.  (17),  the  previously  determined 
flows  per  unit  width  at  time  level  (t+1 ) are  used  to  establish  the  spatial  flow 
gradients.  Thus  in  the  solution  process,  the  computed  flows  and  the  water  levels 
are  always  determined  one-half  time  step  apart,  although  results  are  usually 
presented  only  at  even  time  levels. 

The  difference  formulations  presented  above  for  solution  of  the  three  basic 

unknowns  are  applicable  to  situations  where  grid  elements  can  be  considered  as 

"open  water"  cells  surrounded  by  free  flowing  water  on  all  four  sides.  Usually 

the  majority  of  the  computational  elements  used  to  simulate  a prototype  system 

are  of  this  type,  however  there  are  certain  special  cases  where  the  cross-product 

q 9q  q 3q 

terms  associated  with  convective  acceleration  and  jp  ^ - must  be 

formulated  in  a manner  that  is  consistent  with  the  conditions  that  exist  in  the 

prototype.  For  example,  where  an  impermeable  boundary  such  as  an  island  or  dike 

lies  between  two  adjacent  water  cells  as  illustrated  below  by  the  heavy  line, 

the  flow  gradient  3q  /3y  expressed  as  [qv(i,  j+1)  - qv(i.  j-l)]/2  Ay.  eq.  (6), 

has  little  meaning  since  the  two  flows  used  to  approximate  the  gradients  are  not 

linked  by  a continuous  water  surface. 

Y 


I 


1 
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In  this  case,  a more  meaningful  approximation  of  the  gradient  would  be  [qx(i,j) 
- qx(i,  J-l )]/4y.  Similarly  if  the  impermeable  barrier  were  along  the  bottom 
of  CELL  (i,j),  the  quantity  [qx(i,  j+l)  - qx(i>j)]/Ay  would  define  the  flow 
gradient.  For  water  cells  lying  between  impermeable  barriers,  the  lateral 
flow  gradients  are  set  equal  to  zero.  These  alterations  to  the  finite  differ- 
ence approximations  of  the  cross-product  terms  for  both  the  x and  y directions 
are  made  in  the  model  for  all  possible  combinations  of  impermeable  barrier 
configurations. 

To  facilitate  the  use  of  the  proper  forms  of  the  cross-product  differ- 
ence approximations  in  the  model,  a special  cell  flagging  scheme  has  been 
devised.  Each  element  in  the  grid  is  assigned  a two  digit  number  which 
designates  the  calculations  that  are  to  be  performed  in  describing  the  flow 
gradients  across  the  element  in  each  of  the  two  coordinate  directions.  The 
two  digits  may  each  range  from  one  to  four  with  the  first  digit  applying  to 
the  x-direction  and  the  second  to  the  y-direction.  Under  this  system,  the 
following  finite  difference  approximations  to  the  two  cross-product  terms 
are  made  in  the  model  according  to  the  value  of  the  individual  digits. 


DIGIT  VALUE 

1 

2 

3 

4 


TYPE  OF  APPROXIMATION 

Centered  Difference 
Zero 

Forward  Difference 
Backward  Difference 


For  each  of  the  two  directions,  the  cross-product  terms  are  expressed  as 
follows: 
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Special  Boundary  Conditions 

Because  of  the  complex  geometry  and  flow  behavior  in  most  tidal ly 
influenced  prototype  systems,  there  are  several  different  types  of  special 
flow  boundaries  that  must  be  described  to  properly  simulate  prototype  con- 
ditions with  HYDTID.  For  the  model  application  described  herein,  boundary 
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conditions  can  be  grouped  according  to  the  following  seven  general  classi- 
fications: 

1.  Water-Land.  Boundaries 

2.  Submerged  Barrier  Boundaries 

3.  Overtopping  Barrier  Boundaries 

4.  External  Flow  Boundaries 

5.  Artificial  Offshore  Tidal  Boundaries 

6.  Artificial  Tidal-Storage  Boundaries 

7 . Inundation  Boundaries 

Water-Land  Boundaries  occur  primarily  around  the  perimeter  of  the  system 
although  other  internal  physical  features  such  as  islands,  spoil  banks,  dikes 
and  jetties  can  provide  the  equivalence  of  water-land  boundaries.  A necessary 
requirement  for  computation  of  tidal  hydrodynamics  adjacent  to  these  boundaries 
is  that  the  component  of  flow  normal  to  the  boundary  be  equal  to  zero,  qn  = 0. 
In  the  model,  this  condition  is  handled  by  specifying  the  following  equations 
for  those  elements  where  impermeable  boundaries  exist. 


l*+1(1,j) 

= 0 

(18) 

ly+l(l.j) 

= 0 

(19) 

Submerged  Barrier  Boundaries  are  used  in  the  model  to  describe  interior 
flows  across  such  barriers  as  submerged  reefs,  spoil  banks,  or  pipe  lines  and 
to  provide  flow  control  through  constrictions  such  as  narrow  channels  and 
passes.  An  important  characteristic  of  this  type  of  flow  boundary  is  that 
the  water  levels  on  each  side  of  the  barrier  always  exceed  the  barrier  crest 
elevation.  Under  this  situation,  the  flow  across  the  barrier  can  be  described 
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in  a manner  analogous  to  that  used  for  submerged  weirs,  i.e. 


qn  = ± CB  dB  /g|hi  - h2 | (20) 

where  dg  is  the  depth  of  water  over  the  crest  of  the  weir,  hi  and  h2  are  the 
water  levels  on  each  side  of  the  weir,  and  Cg  is  an  appropriate  discharge 
coefficient. 


Application  of  this  boundary  condition  in  the  model  is  accomplished 
using  the  same  set  of  finite  difference  equations  previously  developed  for 
open  water  cells  except  that  in  place  of  a bottom  resistance  term  involving 
the  energy  slope  in  eq.  (9), a new  term  is  substituted  which  incorporates  the 
effect  of  weir  control.  With  this  substitution  the  only  expressions  affected 
are  those  for  C^x  and  C^.  , eqs.  (11)  and  (13).  For  a typical  barrier  parallel 
to  the  y-axis  and  positioned  at  the  right  side  of  CELL  (i,j),  q^(i,j)  can  be 
determined  using  eq.  (10)  and  the  following, 


C 


fx 


1.0  + 


At  qj_1 (i+1 , j)  - q*_1(i-1.j) 
*7  2"E 


+ 


At  dx  q£_1(i  ,j) 

~ <"CB  dB>2 


in  which 

dB  = 1/2  [ht(i+l , j)  + ht(i .j)]  - zB 

and  zB  is  the  barrier  crest  elevation. 


(21) 
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Similarly  for  a submerged  barrier  parallel  to  the  x-axis  at  the  top  of 
CELL  (i,j). 


1.0 


4t  j*l)  - qJ'’H,J-l) 

1,  2 4y 


+ at  dyl"v~1|1-i)l 

4)1  I',V! 


where  dg  = 1/2  [ht(i,  j+1 ) + ht(i,j)]  - zg 


(22) 


Again,  the  expressions  for  Cfx  and  Cfy  are  essentially  the  same  as  used  by 
Reid  and  Bodine  [1968]  except  that  a term  accounting  for  the  convective 
acceleration  has  been  added.  This  discharge  coefficient  Cg  is  also  the  same 
varying  between  0.01  and  1.3  and  is  determined  in  part  during  the  calibration 
process. 

Overtopping  Barrier  Boundaries  are  included  in  the  model  to  describe 
the  flows  which  occur  periodically  over  dikes,  jetties,  bridge  embankments, 
levees,  etc.  as  water  levels  rise  and  fall  under  the  influence  of  tidal 
action,  winds,  flood  waves,  or  other  external  forces.  Normally  the  water 
level  on  only  one  side  of  these  type  boundaries  exceeds  the  barrier  crest 
elevation  and  under  such  conditions,  flows  are  approximated  in  the  model 
using  an  expression  similar  to  the  broad-crested  weir  formula. 


qn  = ± Cg  AH  i/gAH  (23) 

where  Cg  is  a discharge  coefficient  on  the  order  of  0.5  and  AH  is  the  elevation 
of  the  upstream  total  energy  line  above  the  weir  crest.  For  an  overtopping 
barrier  parallel  to  the  y-axis  and  at  the  right  side  of  CELL  (i,j),  eq.  (23) 
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can  be  written  as, 


q'+,(i.J)  ■ 1 C'dB/gT7  (24) 

"here  dB  ’ hma«  - 2B  " 

and  hm=,v  = ht(i,j)  or  ht(i+l,  j),  whichever  is  greater. 

lilaX 

Similarly  for  a barrier  at  the  top  of  CELL  (i,j), 

qj+1(i,j)  = ± Cg  dB  (25) 

where  dB  = hmax  - ZB  3 

and  hm=i„  = ht(i,j)  or  ht(i,  j+1),  whichever  is  greater. 

As  the  downstream  water  level  rises  due  to  an  increasing  flow  across 
the  barrier,  it  is  possible  that  complete  submergence  of  the  barrier  can 
occur.  When  this  takes  place,  flow  control  for  the  barrier  is  transferred 
in  the  model  to  the  appropriate  submerged  weir  relation  given  by  eq.  (20) 
above.  Therefore  the  overtopping  barrier  is  subject  to  dynamic  boundary 
conditions  according  to  the  water  levels  which  exist  on  both  sides  of  the 
barrier. 

External  Flew  Boundaries  are  specified  in  the  model  where  river  inflows, 
diversions,  and  return  flows  enter  or  leave  the  perimeter  of  the  system. 

These  flows  can  be  expressed  quantitatively  by  the  equation,  qn  = Q(t),  which 
allows  the  external  flows  to  be  introduced  or  withdrawn  in  any  time  dependent 
manner  or  at  a constant  rate  over  a prescribed  period  of  time  at  any  boundary 
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location.  Thus  in  the  model , 


q*+1(i,j)  = Q(t+1 ) (26) 

and 

qj+1(i.j)  = Q(t+1 ) (27) 

where  Q is  the  prescribed  inflow. 

Aside  from  being  used  to  define  inflows  or  withdrawals  along  the  system 
perimeter,  external  flow  boundaries  are  also  necessary  to  describe  flow  con- 
ditions along  the  outer  boundaries  of  fine-grid  sub-models  of  localized  areas. 
This,  of  course,  assumes  that  such  flows  are  known  and  can  be  prescribed  with 
sufficient  detail  in  both  time  and  space  to  properly  define  prototype  conditions. 
When  this  information  is  not  available  from  field  measurements,  a coarse  grid 
model  of  the  entire  system  can  be  used  to  generate  these  internal  flows,  which 
upon  interpolation,  provide  the  necessary  external  inputs  to  the  more  detailed 
fine  grid  model . 

Artificial  Offshore  Tidal  Boundaries  are  used  in  the  model  to  describe 
the  tidal  action  which  occurs  immediately  offshore  of  the  prototype  system. 

These  boundaries  must  be  accurately  defined  since  they  represent  the  boundary 
along  which  the  major  forcing  function  for  excitation  of  the  tidal  hydro- 
dynamics must  be  applied.  In  HYDTID,  tidal  boundaries  are  specified  along  an 
imaginary  line  offshore  and  generally  parallel  to  the  coastline.  The  tidal 
flow  across  this  boundary  corresponding  to  the  excitation  tide  is  computed 
according  to  the  relationship, 

q„  = CT  (Hg  " h>  (28) 
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where  Cj  is  an  admittance  coefficient  taken  approximately  equal  to  the  speed 
of  a gravity  wave  (/g  d) , Hg  is  the  forcing  function  or  known  water  surface 
time  history  specified  at  the  seaward  boundary,  and  h is  the  actual  water 
level  at  the  outer  boundary  of  the  tidal  system.  This  relation  is  consistent 
with  that  suggested  by  Reid  and  Bodine  [1968]. 

For  a typical  computational  element  in  the  model  with  an  ocean  boundary 
on  its  right  side,  the  tidal  flow  can  be  computed  as  follows, 

qj+1(i,j)  = ± CT  [Hg  - ht(i+l , j)]  (29) 

Similarly  for  an  ocean  boundary  at  the  top  of  an  element, 

qy+1(i,j)  = ± CT  [Hg  - ht(i+l , j)]  (30) 

In  these  equations,  the  sign  is  such  that  positive  flow  is  into  the  tidal 

system. 

Artificial  Tidal-Storage  Boundaries  are  specified  in  the  model  to  provide 
the  proper  exchange  of  water  across  external  model  boundaries  which  lie  inside 
the  mouth  of  a tidal  system  where  tidal  flows  occur.  To  eliminate  the  necessity 
to  model  the  entire  water  area  included  in  the  tidal  system,  the  storage  pro- 
vided by  external  water  bodies  adjacent  to  the  modeled  area  can  be  simulated 
by  using  a relation  similar  to  eq.  (28)  applied  along  the  model  perimeter.  In 
this  case,  the  forcing  function,  H^,  is  the  tide  which  occurs  along  the  model 
perimeter  in  the  prototype  system.  Since  this  tide  is  normally  not  available 
from  measured  data,  it  is  often  necessary  to  determine  its  shape  and  phasing 
during  the  verification  of  the  model  by  comparing  computed  and  measured  flows 
inside  the  modeled  area.  Once  these  tidal  conditions  are  established,  cor- 
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relations  can  then  be  developed  with  other  known  tides  which  are  always  avail 
able  from  prototype  records. 


Inundation  Boundaries  are  used  in  the  model  to  account  for  flooding  of 
areas  with  elevations  between  mean  low  and  high  water  levels.  Computations 
at  this  type  of  flow  boundary  are  characterized  by  a check  of  the  water  levels 
in  adjacent  cells  relative  to  the  ground  elevation  to  determine  if  inundation 
is  possible.  If  not,  the  cell  boundary  remains  a Uater-Land  Boundary  and  a 
condition  of  zero  flow  is  specified.  If  flooding  is  possible,  the  boundary 
is  treated  as  if  it  were  an  open  water  cell  and  the  flows  are  computed  according 
to  eqs.  (10)  and  (12). 

Because  of  the  variety  of  different  boundary  conditions  that  must  be 
handled  within  HYDTID , it  is  both  convenient  in  setting  up  the  model  and 
computationally  expedient  to  define  the  different  types  of  cells  a priori. 

Each  cell  in  the  grid  system  is  therefore  identified  or  "flagged"  according 
to  the  flow  conditions  along  its  top  and  right  boundaries.  Figure  3.  (The 
designation  of  the  cell  flags  in  Figure  3 is  completely  independent  of  the 
two-digit  flagging  scheme  used  to  handle  the  cross-product  convective  acceler- 
ation terms. ) 

Although  the  basic  HYDTID  model  described  herein  accounts  for  forty- 
five  combinations  of  different  boundary  conditions  as  defined  along  the  top 
and  right  sides  of  the  cells  only  twenty-eight  were  actually  used  In  modeling 
Masonboro  Inlet.  The  remaining  unused  cells  result  from  permutations  of  the 
seven  basic  boundary  conditions  and  experience  with  applications  of  HYDTID 
in  other  systems.  It  is  probable  that  these  additional  flags  will  find  use 
in  other  modeling  efforts. 


24 


Figure  3 

CELL  FLAGGING  SCHEME  FOR  IDENTI- 
FYING FLOW  BOUNDARY  CONDITIONS 
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Choice  of  Time  Step  and  Mesh  Size 


In  an  explicit  solution  scheme,  the  mesh  size  and  computational  time 
step  are  related  through  mathematical  considerations  arising  from  stability, 
convergence  and  compatibility.  Specifically,  the  following  criterion  must 
be  maintained  for  a stable  solution  to  eqs.  (10),  (12),  and  (17) 


At  < 


AS 


g d, 


(31) 


max 


where  At  is  the  time  step;  As  is  the  mesh  size;  and  d_.„  is  the  maximum  water 
depth  used  in  the  model.  Ideally,  As  (or  Ax  or  Ay)  is  made  as  small  as  possible 
to  provide  maximum  spatial  resolution  and  detail.  This,  however,  can  only  be 
done  at  the  expense  of  computer  storage  and  time.  Thus,  the  selection  of  the 
mesh  size  must  be  based  on  available  prototype  data,  required  model  detail, 
desired  resolution  of  results,  computer  time  and  storage,  as  well  as  the  mathe- 
matical constraint  imposed  by  eq.  (31). 


Implementation  Procedure 

HYDTID  is  designed  so  that  tidal  computations  begin  at  the  seaward 
boundary  and  progress  shoreward  explicitly  over  the  entire  computational 
matrix.  Time  is  advanced  incremently  and  the  computations  are  repeated 
until  the  desired  simulation  period  is  completed.  To  apply  and  use  HYDTID 
the  following  basic  steps  are  required: 


1.  Define  model  resolution  (in  accordance  with  procedures  outlined 
in  the  following  Chapter). 

2.  Determine  physical  boundaries  of  the  area  to  be  modeled  and  super- 
impose a mesh  of  cells  (scaled  to  provide  the  desired  resolution) 
over  the  corresponding  hydrographic  chart  of  the  model  area. 
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3.  Determine  average  MSL  water  depth  (or  land  elevation)  for  each 
cell  in  the  mesh. 

4.  Determine  cell  computational  identification  number  (flag)  according 
to  Figure  3 for  each  element  in  the  model.  Consider  only  the  top 
and  right  boundaries  of  individual  cells  when  establishing  flags. 

5.  Assign  appropriate  two-digit  flag  numbers  to  those  cells  where 
convective  cross-product  terms  must  be  evaluated  by  expressions 
other  than  central  difference  approximations.  (Flags  "11"  need 
not  be  specified.) 

6.  Assign  discharge  coefficients  and  crest  elevations  to  all  sub- 
merged and  overtopping  barriers  in  the  model. 

7.  Estimate  bottom  friction  in  terms  of  Manning  "n"  coefficients 

for  each  water  cell  and  for  each  land  cell  likely  to  be  inundated. 

8.  Based  on  availability  of  tidal  and  hydraulic  data,  select  periods 
for  model  calibration. 

9.  Assemble  all  model  input  data  according  to  the  order  and  formats 
specified  in  the  Appendix. 

10.  Operate  the  model  for  the  selected  calibration  period  until  repe- 
titive conditions  from  one  tidal  cycle  to  the  next  are  achieved 

and  compare  computed  tides  and  velocities  (as  well  as  total  exchange) 
against  measured  data  at  common  sections. 

11.  Adjust  and  refine  the  model  through  bottom  friction,  discharge 
coefficients  and  elevations  of  submerged  barriers  used  to  control 
flow,  and  the  admittance  coefficients  for  excitation  tides  until 
desired  agreement  is  obtained  between  computed  and  measured  results. 

12.  Carry  out  operational  runs  using  the  verified  model.  Depending  on 
the  character  of  any  proposed  changes  in  the  operational  runs,  it 
is  usually  necessary  to  change  some  cell  flag  identifications  to 
account  for  the  corresponding  changes  in  the  bathymetry  or  physical 
features  in  the  system. 
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III.  MODEL  APPLICATION  AND  VERIFICATION 


Prototype  Information 

For  the  12  September  1969  verification  period,  physiographic  and  bathy- 
metric information  at  Masonboro  Inlet  was  obtained  from  a 1"  = 500*  hydrographic 
chart  provided  by  CERC  and  based  on  data  collected  during  September  and  October 
of  1969.  The  chart  contained  two-foot  contours  describing  bottom  topography 
below  MLW  and  spot  elevations  for  ground  levels  above  MHW.  The  datum  for  all 
elevations  was  given  as  MLW.  Also  included  on  the  chart  were  the  jetty  align- 
ment and  the  locations  of  several  recording  tide  gages  and  velocity  stations 
where  data  were  collected  for  the  September  verification  period. 

Seven  tide  gages  were  used  to  record  continuous  tidal  information  in  the 
vicinity  of  Masonboro  Inlet  for  most  of  the  period  0300  to  2100  on  12  September 
1969.  The  spatial  coverage  provided  by  the  gages  is  illustrated  in  Figure  4 
and  is  considered  adequate  for  model  development,  particularly  since  tidal 
levels  at  these  locations  differ  by  only  a few  tenths  of  a foot  and  phase 
changes  are  generally  less  than  an  hour.  For  the  verification  period,  tides 
were  nearly  sinusoidal  in  shape  with  a semi-diurnal  period  on  the  order  of 
12.5  hours.  Representative  of  the  tides  measured  at  all  of  the  gages  is  that 
for  Gage  0 located  on  the  ocean  side  oi  Wrightsville  Beach  north  of  the  jetty. 
This  offshore  tide,  shown  in  Figure  5,  provides  the  primary  tidal  influence  to 
the  Masonboro  Inlet  area,  and  as  such,  is  of  major  importance  for  providing 
tidal  excitation  in  the  model. 


29 


Also  shown  in  Figure  4 are  five  ranges  where  current  velocity  measure- 
ments were  made  during  the  verification  period.  For  each  range,  the  velocity 
measurements  generably  covered  the  period  0700  to  2100  on  12  September  1969. 
Three  stations  spaced  from  150  to  300  feet  apart  were  used  to  provide  a 
description  of  velocity  variations  across  each  range.  Also  for  all  of  the 
stations,  plots  of  velocity  vs.  time  at  the  surface,  mid-depth,  and  bottom 
were  provided  for  most  of  the  verification  period. 

Information  describing  the  weir  section  of  the  jetty  was  obtained  from 
CERC  and  used  in  structuring  the  model  of  the  inlet  area.  The  weir  section 
of  the  jetty  extends  approximately  1000  feet  seaward  from  the  north  back  of 
the  inlet,  and  has  a crest  elevation  of  2.0  feet  above  MLW,  a height  that  is 
exceeded  during  periods  of  normal  high  tide.  Beyond  the  weir  section,  the 
jetty  continues  approximately  1900  feet  offshore.  For  purposes  of  this 
study,  it  has  been  assumed  that  this  section  of  the  jetty  is  never  overtopped 
by  normal  tides  and  is  completely  impermeable. 

Model  Resolution 

As  noted  previously,  selection  of  the  mesh  size  for  application  of  HYDTID 
to  Masonboro  Inlet  is  a function  of  the  system  geometry,  the  spatial  density 
of  prototype  data  measurements,  the  desired  resolution  of  the  model  results, 
and  the  economics  of  the  model  operation.  Of  course,  the  storage  and  compu- 
tational efficiency  of  available  computer  hardware  must  also  be  considered. 

The  final  selection  of  a mesh  size  in  this  study  has  been  based  on  an  evaluation 
of  each  of  these  factors  for  the  Masonboro  Inlet  problem,  as  well  as,  on  experi- 
ence in  applying  HYDTID  to  several  different  tidal  systems. 
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A review  of  the  Masonboro  Inlet  study  area.  Figures  1 and  4,  including 
the  configuration  of  islands,  inter-connected  channels,  and  inundated  areas, 
suggests  that  a mesh  size  on  the  order  of  1000  to  1200  feet  would  provide 
sufficient  resolution  for  generally  describing  these  features  and  allowing 
some  simulation  of  tidal  flows  and  water  surface  fluctuations.  While  this 
level  of  resolution  would  not  permit  the  detailed  description  of  flows  and 
velocities  within  channels  themselves,  it  would  provide  for  reasonably  accurate 
simulation  of  the  gross  flow  distribution  at  the  major  channel  junctures  and 
through  the  main  inlet. 

However  to  be  consistent  with  the  spatial  resolution  of  available  proto- 
type data,  and  in  particular  the  velocity  measurements,  a mesh  size  of  100  to 
300  feet  more  nearly  conforms  to  the  spacing  between  individual  velocity 
stations.  Additionally,  this  resolution  would  allow  a better  representation 
of  the  variable  bottom  topography  in  the  inlet  area,  thereby  permitting  more 
accurate  simulation  of  localized  flows  within  the  major  branches  of  the  system. 

The  advantages  of  using  still  finer  resolution  in  the  model  of  Masonboro 
Inlet  are  obvious,  however  as  is  often  the  case,  it  is  the  cost  of  operating 
the  model  that  ultimately  determines  the  fineness  with  which  physical  detail 
can  be  described.  A reductior  in  the  mesh  size  results  not  only  in  more 
computational  elements  to  cover  a iven  area,  but  also  requires  a reduction 
in  the  computational  time  step  because  of  the  stability  criterion,  eq.  (31), 
imposed  on  the  solution  scheme.  For  example  if  the  element  size  is  reduced  by 
half,  the  number  of  calculations  required  to  simulate  a specified  period  of 
real  time  is  theoretically  increased  by  a factor  of  eight,  thus  increasing  the 
cost  of  operating  the  model  also  by  a factor  of  about  eight. 
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To  aid  in  the  selection  of  a near-optimal  mesh  size  for  the  Masonboro 
Inlet  model,  the  computer  time  required  to  operate  models  with  different 
element  sizes  ranging  from  200  to  1200  feet  was  estimated  for  simulation 
periods  of  one  tidal  cycle  (12.5  hours).  The  study  area  identified  in  Figure 
1 was  used  for  structuring  these  hypothetical  models,  and  experience  with 
HYDTID  in  other  systems  was  used  to  estimate  the  computer  time  requirements 
for  each  of  the  selected  mesh  sizes. 

Figure  6 summarizes  the  UNIVAC  1108  CPU  time  requirements  for  HYDTID 
models  of  several  different  tidal  systems  with  which  WRE  has  had  modeling 
experience.  The  abscissa  in  Figure  6 is  equal  to  the  product  of  the  number 
of  water  cells  included  in  the  model  where  calculations  must  be  performed  and 
the  number  of  time  steps  required  to  simulate  24  hours  of  real  time.  For  ex- 
ample, if  a mod°l  consisted  of  100  water  cells  and  used  a time  step  of  one 
minute,  a total  of  (24  x 60)/l  = 1440  time  steps  would  be  required  to  simu- 
late one  day,  and  the  abscissa  value  used  in  Figure  6 would  be  equal  to 
100  x 1440  = 144000.  Thus,  knowing  the  size  of  the  prototype  system,  including 
its  maximum  depth,  and  selecting  a trial  mesh  size,  the  allowable  time  step  can 
be  determined  from  eq.  (31),  and  finally,  the  computer  time  and  cost  require- 
ments for  operating  the  model  can  be  estimated  from  Figure  6.  The  magnitude 
of  the  error  which  might  be  expected  in  such  a time  estimate  is  indicated  by 
the  scatter  of  the  points.  Since  not  all  tidal  systems  are  exactly  alike, 
model  operation  times  also  vary,  not  only  because  of  mesh  size  and  time  step 
differences,  but  also  because  of  the  number  and  variety  of  boundary  conditions 
required  by  models  of  different  systems.  It  should  also  be  noted  that  the 
UNIVAC  1108  CPU  times  used  in  developing  Figure  6 were  all  determined  based  on 
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UNIVAC  1108  CPU  TIME  FOR  24  HOUR  SIMULATION,  MINUTES 


Figure  6 

COMPUTER  TIME  REQUIREMENTS  FOR  'HYDTID'  MODEL 
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the  Tracor,  Inc.  computer  system  in  Austin,  Texas,  and  results  from  this  figure 
may  not  be  completely  consistent  with  other  UNIVAC  1108  systems. 

For  the  Masonboro  Inlet  study  area,  the  results  of  applying  Figure  6 with 
different  mesh  sizes  are  summarized  in  Table  1. 


Table  1.  Model  Characteristics  for  Masonboro  Inlet 


Mesh 

Size 

(feet) 

Grid  Size 

Total 

No.  of 
Cells 

Total  No. 
of  Water 
Cells 

Time 

Step 

(Sec.) 

Computer  Time 
per  Tidal 
Cycle  (min.) 

200 

72  x 96 

6912 

4896 

4.2 

544.0 

300 

48  x 64 

3072 

2176 

6.3 

161.6 

400 

36  x 48 

1728 

1224 

8.3 

68.0 

600 

24  x 32 

768 

544 

12.5 

20.2 

800 

18  x 24 

432 

306 

16.6 

8.5 

1200 

12  x 16 

192 

136 

25.0 

2.5 

This  table  illustrates  that  to  model  Masonboro  Inlet  using  an  element  size 
consistent  with  the  spacing  of  velocity  stations  (100-300  feet)  requires  con- 
siderable computer  time.  Even  if  mesh  sizes  on  the  order  of  400  to  600  feet 
were  used  to  structure  the  model  , the  advantage  of  these  models  in  terms  of 
computer  cost  savings  would  likely  be  outweighed  by  their  inability  to  adequately 
simulate  localized  flows  with  the  coarser  resolution.  Thus  it  does  not  appear 
feasible  to  model  the  entire  Masonboro  Inlet  area  with  a single  grid  configuration 
If  the  mesh  size  must  be  consistent  with  the  spatial  density  of  available  proto- 
type data. 

An  alternate  procedure  is  to  describe  the  hydrodynamics  of  the  entire  area 
on  a large  scale  using  a coarse  grid  model,  and  then  use  the  simulated  results 
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from  this  model  to  provide  the  external  inputs  to  a fine  grid  sub-model  of  a 
localized  area  where  more  detail  is  desired.  This  dual-model  approach  results 
in  a computer  cost  savings  if  the  combined  operation  time  of  the  two  models  is 
less  than  that  of  a single  fine  grid  model  of  the  entire  area.  For  the  Masonboro 
Inlet  area,  this  can  be  shown  to  be  the  case  by  considering  a coarse  grid  model 
of  the  entire  area  based  on  an  element  size  of  1200  feet  and  a fine  grid  sub- 
model of  only  the  inlet  area  using  a mesh  size  of  300  feet.  The  extent  of 
these  models  is  shown  in  Figure  7 and  covers  the  same  area  as  that  designated 
in  Figure  1 as  the  study  area. 

From  Table  1,  the  computer  time  required  to  simulate  one  tidal  cycle 
with  the  1200'  coarse  grid  model  should  be  about  two  to  three  minutes.  Using 
a similar  procedure,  computer  time  for  the  fine  grid  sub-model  can  be  estimated 
from  Figure  6 as  follows: 

Mesh  Total  Total  No.  Time  Computer  Time 

Size  Grid  Size  No.  of  of  Water  Step  per  Tidal 

(feet)  Cells  Cells  (Sec.)  Cycle  (min.) 

300  33  x 30  990  360  5.0  33.6 

Hence  the  combined  computer  time  requirement  of  the  two  models  is  on  the  order 
of  36  minutes  and  less  than  the  162  minutes  given  in  Table  1 for  a single  model 
covering  the  entire  study  area  with  a 300'  mesh  size.  Based  on  these  projections, 
the  model  of  Masonboro  Inlet  has  been  developed  using  the  dual -model  approach 
with  mesh  sizes  of  1200  and  300  feet. 
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Dual-Model  Procedure 


Implementation  of  the  dual-model  concept  for  verification  of  the 
12  September  1969  conditions  in  Masonboro  Inlet  involved  first  the  cali- 
bration of  the  coarse  grid  model  for  generalized  flow  simulation  and 
reproduction  of  tidal  amplitudes,  and  secondly,  the  calibration  of  the 
fine  grid  sub-model  to  describe  localized  flows  and  velocities.  While 
the  entire  verification  process  required  the  development  and  tuning  of 
two  hydrodynamic  models,  calibration  of  the  coarse  grid  model  resulted 
in  information  directly  usable  in  the  fine  grid  sub-model.  Linkage 
between  the  two  models  was  provided  by  the  flows  computed  in  the  coarse 
grid  model  across  those  cell  boundaries  which  coincided  with  the  external 
limits  of  the  fine  grid  sub-model.  Initial  water  levels  in  the  fine  grid 
sub-model  were  determined  by  interpolation  of  corresponding  water  levels 
computed  with  the  coarse  grid  model. 

An  important  consideration  in  implementing  the  dual-model  concept 
Is  that  the  relation  between  depth  and  volume  be  the  same  in  the  two 
models  for  that  portion  of  the  prototype  contained  in  both.  Specifically, 
this  means  that  ground  elevations  between  MLW  and  MHW  must  be  assigned 
to  the  individual  cells  in  each  of  the  models  in  a consistent  manner  so 
that  for  any  water  level  above  MLW,  the  amount  of  inundation  is  essenti- 
ally the  same, with  both  models  providing  equal  quantities  of  water  storage 
capacity.  If  the  depth-volume  relations  are  not  the  same,  then  computed 
tidal  levels  in  the  models  will  not  agree  when  boundary  flows  to  the  two 
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systems  are  equal.  The  degree  to  which  the  depth-volume  relations  in  the 
two  models  can  be  made  to  coincide  depends  on  the  relative  difference 
between  the  grid  mesh  sizes  and  the  variation  of  the  topography  between 
the  MLW  and  MHW  levels. 

It  should  be  pointed  out  that  the  differences  in  the  computed  hydro- 
dynamics which  result  from  not  having  compatible  bottom  depths  and  ground 
elevations  in  the  two  models  are  most  discernible  in  the  tidal  levels 
rather  than  in  the  channel  flows  and  velocities.  The  main  reason  for 
this  is  that  simulated  flows  from  the  coarse  grid  model  are  used  to  excite 
the  fine  grid  sub-model,  and  hence,  flows  within  the  fine  grid  system 
directly  reflect  these  inputs.  If  the  differences  which  occur  in  com- 
puted tidal  levels  due  to  inconsistent  depth-volume  relations  are  small 
compared  to  the  total  depths,  the  channel  cross-sectional  areas  will  not 
change  significantly  and  consequently,  the  velocities  will  be  only  slightly 
affected . 

To  minimize  the  effects  of  inconsistent  depth-volume  relations  in 
the  two  models,  the  bottom  and  ground  elevations  should  first  be  estab- 
lished for  each  cell  in  the  fine  grid  sub-model,  and  then, based  on  this 
Information  together  with  the  inlet  geometry,  appropriate  average  ele- 
vations can  be  assigned  to  each  cell  in  the  coarse  grid  model.  While 
the  order  of  this  procedure  Is  the  reverse  of  that  normally  used  In 
model  calibration  and  verification,  it  provides  for  better  compati- 
bility between  the  depth-volume  relations  even  though  a few  small 
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depth  adjustments  may  be  necessary  during  subsequent  calibration  of  the 
models.  Aside  from  this  variation,  the  development  of  the  two  models 
followed  the  general  implementation  procedure  outlined  in  the  previous 
chapter. 

Coarse  Grid  Model.  The  coarse  grid  configuration  used  to  describe 
the  physiography  of  the  Masonboro  Inlet  area  is  shown  in  Figure  8.  Using 
a 1200'  grid  overlay  with  the  1"  = 500'  hydrographic  chart,  water- land 
boundaries  were  established  to  generally  conform  to  the  MHW  line.  In- 
cluded in  the  grid  network  are  135  water  cells  covering  an  area  of  about 
seven  square  miles.  The  shaded  cells  in  Figure  8 either  represent  dry 
land  where  inundation  never  occurs  under  normal  tidal  conditions  or  they 
denote  "dead"  cells  where  no  calculations  are  required  but  are  needed  in 
the  model  to  acconmodate  the  two-dimensional  Indexing  scheme  and  to  facili- 
tate the  computations  in  certain  boundary  cells.  It  should  be  noted  that 
an  extra  row  and  column  of  "dead"  cells  are  Included  in  the  grid  along 
the  top  and  right  sides,  respectively.  These  additional  cells  permit 
definition  of  certain  key  variables  that  are  used  in  the  calculation 
process. 

Bottom  and  ground  elevations  were  assigned  in  the  coarse  grid  model 
using  the  general  procedure  outlined  above  in  an  effort  to  minimize  the 
differences  In  the  depth-volume  relations  of  the  two  models.  Even  using 
this  procedure,  inconsistencies  reflected  In  the  subsequent  model  results 
could  not  be  avoided  because  of  the  coarse  grid  resolution.  Using  the 
1"  * 500'  hydrographic  chart,  MLW  elevations  were  first  assigned  to  the 
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cells  in  the  fine  grid  sub-model,  and  then  appropriate  averages  were  deter- 
mined for  the  coarse  grid  cells.  For  coarse  grid  cells  outside  the  external 
limits  of  the  fine  grid,  MLW  elevations  were  assigned  directly  from  the  hydro- 
grapnic  chart  using  a 1200'  grid  overlay.  These  MLW  elevations  were  then 
converted  to  a MSL  datum  by  applying  a uniform  vertical  shift  of  1.3  feet, 
the  difference  between  MSL  and  MLW  as  reported  by  Tide  Tables  [1969]  for  the 
general  coastal  area  in  the  vicinity  of  Masonboro  Inlet.  The  field  of  MSL 
elevations  used  in  the  coarse  grid  model  to  simulate  the  12  September  1969 
conditions  is  given  in  Figure  9.  Asterisks  in  Figure  9 denote  land  cells  which 
are  never  inundated,  and  have  a MSL  elevation  of  +99  feet  assigned  to  them  in 
the  model . 

To  adequately  describe  all  of  the  physiographic  features  and  flow  con- 
trols included  in  the  Masonboro  Inlet  area,  several  of  the  special  boundary 
conditions  discussed  in  the  previous  chapter  were  used.  Figure  10  shows  the 
locations  of  these  boundaries  at  their  respective  cells  in  the  coarse  grid. 

Cell  boundaries  with  no  designation  in  Figure  10  are  subject  to  the  normal  flow 
of  water  as  described  by  the  basic  equations  of  motion,  eqs.  (1),  (2),  and  (3). 

The  submerged  and  overtopping  barrier  boundaries  identified  in  Figure  10 
were  used  to  provide  flow  control  in  the  major  channels  and  to  describe  flow 
constrictions  which  exist  in  the  prototype  such  as  the  weir  section  of  the 
jetty.  Initial  selection  of  barrier  heights  and  discharge  coefficients  was 
based  on  physiographic  detail  in  the  vicinity  of  the  barriers  and  experience 
using  barriers  in  models  of  other  tidal  systems.  Adjustments  to  these  values 
were  made  during  calibration  of  the  model  to  provide  better  simulation  of  the 
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MSL  BOTTOM  AND  GROUND  ELEVATIONS  USED 
IN  VERIFICATION  OF  COARSE  GRID  MODEL 
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flow  distribution  in  the  system.  The  barrier  heights  and  discharge  coefficients 
finally  adopted  are  given  in  Table  2. 


Table  2. 

Characteristics 

of  Barriers 

Used  in  Coarse  Grid 

Model 

Barrier 

Number 

Cell 

I 

Indices 

J 

Top/Side 

Boundary 

Discharge 

Coefficient 

MSL 

Elevation 

1 

6 

6 

Side 

0.5 

-6.3 

2 

6 

7 

Side 

0.4 

-1 .3 

3 

6 

8 

Side 

0.4 

-2.3 

4 

7 

8 

Side 

0.5 

+0.7 

5 

7 

8 

Top 

1.3 

-12.2 

6 

3 

9 

Side 

0.9 

-4.5 

7 

5 

9 

Side 

0.9 

-5.7 

8 

6 

10 

Top 

0.9 

-8.2 

9 

7 

10 

Side 

0.9 

-7.2 

10 

9 

10 

Side 

0.9 

-10.2 

11 

11 

10 

Top 

0.9 

-3.2 

Tidal  excitation  was  provided  around  the  perimeter  of  the  coarse  grid  model 
using  artificial  offshore  tidal  boundaries  on  the  seaward  side  of  the  inlet  and 
artificial  tidal -storage  boundaries  around  the  inland  periphery  of  the  model. 
Although  designated  separately  in  Figure  10,  these  two  types  of  boundary  con- 
ditions are  handled  identically  in  the  model  with  different  input  tides  specified 
at  various  locations  to  reflect  local  tidal  conditions  throughout  the  inlet  area. 
In  the  coarse  grid  model,  four  different  tides  were  used  to  excite  the  system  and 
assigned  to  appropriate  cells  according  to  Table  3. 

The  four  tides  referred  to  above  are  plotted  in  Figure  11  for  the  veri- 
fication period  0400  to  2100  on  12  September  1969.  Tide  No.  1 corresponds  to 
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Table  3.  Assignment  of  Exciting  Tides  in  Coarse  Grid  Model 


Tidal  Cell 
Number 

Cell 

I 

Indices 

J 

Top/Side 

Boundary 

Tidal  Identi- 
fication No. 

1 

3 

1 

Top 

1 

2 

4 

1 

Top 

1 

3 

5 

1 

Top 

1 

4 

6 

1 

Top 

1 

5 

7 

1 

Top 

1 

6 

8 

1 

Top 

1 

7 

9 

1 

Top 

1 

8 

10 

1 

Top 

1 

9 

1 

2 

Side 

1 

10 

11 

2. 

Side 

1 

11 

1 

3 

Side 

1 

12 

11 

3 

Side 

1 

13 

1 

4 

Side 

1 

14 

11 

4 

Side 

1 

15 

1 

5 

Side 

1 

16 

11 

5 

Side 

1 

17 

1 

6 

Side 

1 

18 

11 

6 

Side 

1 

19 

1 

7 

Side 

1 

20 

11 

7 

Side 

1 

21 

1 

8 

Side 

2 

22 

11 

8 

Side 

1 

23 

1 

9 

Side 

2 

24 

11 

9 

Side 

1 

25 

1 

10 

Side 

4 

26 

11 

10 

Side 

3 

27 

1 

11 

Side 

4 

28 

11 

11 

Side 

3 

29 

1 

12 

Side 

4 

30 

11 

12 

Side 

3 

31 

11 

13 

Side 

3 

32 

11 

15 

Side 

3 

47 


48 


the  prototype  tide  measured  at  Gage  0 on  the  open  coast,  and  it  is  specified 
in  the  model  at  each  of  the  cells  designated  in  Figure  10  as  artificial  off- 
shore tidal  boundaries.  The  other  three  tides  were  developed  as  part  of 
the  calibration  process,  and  they  represent  the  inputs  necessary  to  simu- 
late the  general  character  of  the  flows  in  the  major  channels  during  the 
verification  period.  In  the  initial  stages  of  model  calibration,  the  Gage 
5 tide  was  specified  as  Tide  No.  2,  the  Gage  4 tide  was  used  for  Tide  No. 

3,  ana  Tide  No.  4 was  set  equal  to  the  Gage  6 tide.  During  calibration 
of  the  model,  these  three  interior  input  tides  were  adjusted  in  both  phase 
and  amplitude  until  the  necessary  flow  simulation  was  obtained  for  input 
to  the  fine  grid  model. 

Recognizing  that  future  applications  of  the  model  to  other  conditions 
would  require  that  all  four  of  these  input  tides  be  known,  a procedure  has 
been  devised  whereby  approximations  of  the  three  interior  tides  (Nos.  2,  3, 
and  4)  can  be  obtained  from  the  offshore  tide  at  Gage  0 (Tide  No.  1).  In 
this  analysis,  it  has  been  assumed  that  Gage  0 is  a permanent  installation, 
and  that  historical  tidal  records  can  always  be  obtained  from  this  station. 
Using  the  calibrated  input  tides  necessary  for  simulation  of  the  12  September 
1969  conditions  (Figure  11),  correlations  have  been  developed  between  each 
of  the  three  interior  tides  and  the  Gage  0 tide.  Because  of  the  similarity 
between  the  tides,  these  correlations  have  been  based  on  applying  appropriate 
linear  scaling  factors  to  the  Gage  0 tidal  amplitudes  lagged  according  to  the 
phase  differences  determined  in  model  calibration.  To  describe  these  tides, 
a bi-modal  Fourier  approximation  of  the  following  form  has  been  used. 
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Si  {Ao  + C o i Sin  [ (2tt/T ) • (t  + At^. ) + 4>0i] 
+ C o 2 Sin  [ ( 4tt/T ) • ( t + Ati)  + $02]}  + HSi 


(32) 


In  this  expression,  the  subscript  i refers  to  the  identification  number  of 
the  input  tide  used  in  the  coarse  grid  model  (1  £ i £ 4),  Hi  is  the  MLW  tidal 
elevation,  Aq  is  the  Gage  0 mean  tide  level  with  respect  to  MLW  averaged  over 
one  tidal  cycle  having  a duration  of  T hours,  t is  the  time  in  hours  after 
the  beginning  of  the  tidal  cycle  when  Hi  occurs,  C0i  and  C02  are  Fourier 
coefficients  determined  for  approximation  of  the  tide  at  Gage  0,  4>0i  and  $02 
are  the  corresponding  Fourier  phase  angles  for  the  Gage  0 tide,  Ati  is  the 
phase  difference  between  the  i^1  tide  and  the  Gage  0 tide.  Si  is  an  appropriate 
linear  scaling  factor  for  the  i**1  tide,  and  HSi  is  the  corresponding  vertical 
tidal  adjustment  for  the  ith  tide. 

For  the  period  0400  to  1650  on  12  September  1969,  the  following  Fourier 
constants  have  been  determined  for  the  Gage  0 tide. 

A = 2.38920 

0 

Co  1 = 2.08673 

4>oi  = -0.10846 

C02  = 0.04220 

*02  = 1.31290 

To  determine  hourly  values  of  the  four  input  tides  used  for  final  verification 
of  the  coarse  grid  model,  the  constants  given  in  Table  4 were  applied  to 
eq.  (32). 
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Table  4.  Tidal  Constants  Used  in  Verification 


of  Coarse  Grid  Model 


Tidal  Identi- 
fication No. 

Si 

HSi 

Ati 

1 

IBM 

|B 

2 

m 

-1.70 

3 

■Pi 

-1.05 

4 

-1 .45 

It  should  be  pointed  out  that  the  general  use  of  input  tides  determined 
for  other  simulation  periods  by  eq.  (32)  and  the  constants  in  Table  4 are 
limited  somewhat  by  the  degree  to  which  these  tidal  conditions  agree  with 
those  of  the  12  September  1969  verification  period.  The  constants  in  Table  4 
were  determined  through  analysis  of  only  one  tidal  period  using  the  tidal 
records  which  were  readily  available.  Hence  a certain  amount  of  caution 
should  be  exercised  when  applying  them  to  other  periods  to  assure  that  at 
least  general  similarity  exists  between  the  new  tides  and  those  of  12  Sep- 
tember 1969.  Should  additional  prototype  data  be  available  for  Masonboro 
Inlet  during  other  periods,  it  would  be  advantageous  to  test  these  constants 
through  operation  of  the  model,  and  possibly  develop  other  correlations  more 
applicable  to  different  tidal  conditions.  For  purposes  of  this  study,  it  ’s 
believed  that  these  constants  provide  adequate  descriptions  of  the  tidal 
inputs  because  of  the  similarities  in  the  actual  measured  tides. 

The  specification  of  boundary  conditions  for  each  cell  in  the  model 
was  handled  using  the  cell  flagging  scheme  described  in  Figure  3.  Using 
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these  cell  Identification  numbers  and  the  boundary  conditions  identified  in 
Figure  10,  flags  were  assigned  to  each  cell  in  the  coarse  grid  model  as  illus- 
trated in  Figure  12.  In  a similar  manner  those  cells  which  required  special 
finite  difference  formulations  of  the  convective  acceleration  cross-product 
terms  were  identified  using  the  two-digit  flagging  scheme  described  earlier. 
Only  those  cells  which  required  flags  other  than  "11"  were  identified  and 
input  into  the  model  as  data  since  the  initialization  of  program  variables  in 
the  model  includes  assigning  "11"  to  all  cells.  The  field  of  two-digit  con- 
vective acceleration  cell  flags  is  shown  in  Figure  13. 

Manning's  "n"  values  were  assigned  to  each  cell  in  the  coarse  grid 
model  and  input  into  the  program  as  part  of  the  basic  cell  data.  Although 
slight  adjustments  were  made  to  the  Manning's  "n"  field  during  calibration 
of  the  model,  initial  values  were  determined  using  the  following  relation- 
ships which  have  been  developed  through  experience  with  modeling  other  tidal 
systems . 

n = 0.0377  + 0.00667  z,  z < -4 

n = 0.0550  + 0.00500  z,  -4  < z < 0 

n = 0.0550  + 0.00500  z,  z > 0 

In  these  relations,  z is  the  MSL  bottom  or  ground  elevation  measured  nega- 
tively downwards  as  defined  in  the  basic  equations  of  motion.  Values  obtained 
from  these  relations  range  from  0.0177  for  water  depths  of  30  feet  to  0.0500 
for  water  depths  of  one  foot. 

Calibration  of  the  coarse  grid  model  involved  the  successive  operation 
of  the  model  under  the  prescribed  input  conditions  of  12  September  1969, 
including  an  average  northeast  wind  with  a velocity  of  four  knots,  and 
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the  appropriate  adjustment  of  input  tides,  barrier  elevations  and  coefficients, 
bottom  and  ground  elevations,  bottom  friction,  and  tidal  discharge  coefficients 
until  simulated  results  agreed  satisfactorily  with  measured  prototype  tide  and 
flow  conditions.  In  accordance  with  the  fundamental  concept  of  applying  the 
dual-model  approach,  major  emphasis  during  calibration  of  the  coarse  grid  model 
was  placed  on  tidal  verification  and  on  simulation  of  the  generalized  flow 
distribution  between  the  major  channels  in  the  Masonboro  Inlet  area.  Velocity 
verification  was  a consideration  although  comparisons  between  simulated  results 
and  measured  prototype  data  were  limited  only  to  the  average  velocity  condition 
across  each  of  the  ranges  since  the  coarse  model  resolution  did  not  provide  for 
direct  comparisons  of  actual  field  measurements. 

Flow  comparisons  were  made  at  each  of  the  interior  velocity  ranges 
identified  in  Figure  4.  Using  channel  cross-sections  determined  from  the 
1"  = 500'  hydrographic  chart,  maximum  ebb  and  flood  flows  in  the  prototype 
were  computed  based  on  the  detailed  velocity  measurements  at  the  three  stations 
across  each  range.  The  plotted  cross-sections  used  are  shown  in  Figure  14. 

From  the  model,  flows  across  the  appropriate  boundaries  of  the  cells  given  in 
Table  5 were  used  for  comparison.  Direct  comparison  of  computed  and  measured 
flows  across  Range  1 was  not  possible  since  a confined  channel  does  not  exist 
at  this  location  and  therefore,  compatible  cross-sectional  areas  could  not  be 
determined  for  both  the  model  and  the  prototype.  Based  on  the  good  agreement, 
however,  between  simulated  flows  and  those  determined  from  prototype  measure- 
ments at  the  other  four  ranges,  it  is  expected  that  the  computed  flows  across 
Range  1 are  also  reasonable.  Figure  15  presents  the  flow  comparisons  at  these 
other  ranges  for  both  the  ebb  and  flood  conditions  which  occurred  at  approxi- 
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Range  2 Range  3 Range  4 Range  5 

Figure  15 

COMPARISON  OF  PROTOTYPE  AND  SIMULATED  FLOWS  FROM  COARSE  GRID  MODEL 


Table  5.  Coarse  Grid  Cells  Where 
Flows  Were  Compared 


Range 

Number 

Cell 

I 

Indices 

J 

Top/Side 

Boundary 

2 

7 

8 

Top 

3 

6 

10 

Top 

4 

9 

10 

Side 

5 

3 

9 

Side 

mately  1100  and  1800  hours,  respectively,  on  12  September  1969.  Based  prima- 
rily on  these  flow  simulations,  the  coarse  grid  model  was  considered  calibrated 
and  sufficiently  verified  to  be  used  to  provide  the  input  flows  to  the  fine 
grid  sub-model. 


Verification  of  computed  tidal  levels  from  the  coarse  grid  model  was 
determined  by  comparison  of  simulated  tides  with  those  measured  at  the  six 
gages  in  the  Masonboro  Inlet  area  identified  in  Figure  4.  The  period  0400  to 
2100  on  12  September  1969  was  used  for  tidal  verification,  and  tides  computed 
at  the  following  cells  in  the  model  were  used  for  comparison. 

Table  6.  Coarse  Grid  Cells  Where 
Tides  Were  Compared 

Tide  Gage  Cell  Indices 

Number  I J 


7 6 

7 8 

6 10 

9 10 

3 9 

7 13 
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1 

2 

3 

4 

5 

6 


Selection  of  these  cell  locations  for  comparison  of  tides  facilitated  not  only 
tidal  verification  but  also  flow  verification  since  computed  results  for 
several  of  these  cells  were  used  in  both  verification  processes.  Because 
only  a limited  number  of  cells  could  be  identified  in  the  model  for  output 
of  hourly  hydrodynamics,  some  of  the  above  cells  may  not  fall  directly  on  tide 
gage  locations,  however  the  error  in  using  these  slightly  shifted  tides  is  not 
appreciable.  Tidal  simulations  using  results  from  the  calibrated  coarse  grid 
model  are  presented  in  Figures  16  through  21  for  the  six  tide  gages,  and  they 
are  considered  adequate  for  purposes  of  this  study. 

Fine  Grid  Sub-Model . The  areal  extent  of  the  fine  grid  sub-model  covers 
the  immediate  inlet  area  and  is  generally  bounded  on  the  inland  side  by  the 
three  interior  velocity  ranges  (3,  4,  and  5)  and  on  the  seaward  side  by  the 
jetty  and  a portion  of  the  ocean.  Figure  7 shows  the  relation  of  the  fine  grid 
area  to  the  coarse  grid  area.  Again  using  the  1"  = 500'  hydrographic  chart  and 
a 300'  grid  overlay,  water-land  boundaries  were  established  in  the  fine  grid 
model,  and  the  computational  grid  network  was  set  up  according  to  that  shown 
in  Figure  22.  Included  in  the  model  are  approximately  360  water  cells  cover- 
ing slightly  over  one  square  mile  of  surface  area.  The  field  of  MSL  bottom 
and  ground  elevations  assigned  to  each  cell  in  the  grid  is  presented  in 
Figure  23.  Again,  asterisks  denote  "dead"  cells  where  computations  are  not 
required  and  where  ground  elevations  are  set  equal  to  +99  feet. 

The  special  boundary  conditions  used  in  the  fine  grid  sub-model  are 
shown  in  Figure  24.  To  be  compatible  with  the  coarse  grid  model,  several 
submerged  barrier  boundaries  were  specified  in  the  fine  grid  sub-model  at 
the  locations  shown.  Barrier  elevations  and  discharge  coefficients  were 
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CONFIGURATION  USED  TO  DESCRIBE  SEPTEMBER  1969  CONDITIONS 


IN  VERIFICATION  OF  FINE  GRID  SUB-MODEL 


initially  prescribed  in  a manner  consistent  with  those  used  in  the  coarse 
grid  model,  and  then  adjusted  during  model  calibration.  Table  7 summarizes 
the  barrier  characteristics. 


Table  7.  Characteristics  of  Barriers  Used 
in  Fine  Grid  Sub-Model 


Barrier 

Number 

Cell 

I 

Indices 

J 

Top/Side 

Boundary 

Discharge 

Coefficient 

MSL 

Elevation 

1 

19 

10 

Top 

-21.2 

2 

20 

10 

Top 

C-90 

-16.3 

3 

17 

17 

Top 

1.00 

-21 .2 

4 

18 

17 

Top 

1.00 

-16.2 

5 

19 

17 

Top 

0.90 

-5.2 

6 

5 

18 

Side 

0.90 

-6.2 

7 

5 

19 

Side 

0.90 

-4.2 

8 

29 

22 

Side 

0.90 

-19.2 

9 

29 

23 

Side 

0.90 

-16.2 

10 

16 

24 

Top 

0.90 

-10.2 

11 

17 

24 

Top 

0.90 

-15.2 

As  stated  earlier,  tidal  flows  determined  from  the  coarse  grid  model  are 
used  to  excite  the  fine  grid  sub-model.  These  flows  are  prescribed  around  the 
perimeter  of  the  fine  grid  sub-model  at  appropriate  locations  where  the  external 
boundaries  of  the  model  intersect  the  various  channels  and  passes  in  the  inlet 
area.  These  locations  are  identified  in  Figure  24  as  external  inflow  boundaries. 

The  determination  of  the  input  flows  for  each  of  these  external  inflow 
boundaries  in  the  fine  grid  sub-model  requires  both  temporal  and  spatial  adjust- 
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merits  to  the  coarse  grid  flows.  First  of  all,  since  the  computational  time 
steps  used  in  the  two  models  are  not  the  same  (At  = 20  seconds  in  the  coarse 
grid  model  and  At  = 5 seconds  in  the  fine  grid  sub-model),  the  coarse  grid 
flows  must  be  interpolated  to  obtain  the  intermediate  values  required  by  the 
fine  grid  sub-model.  Also  because  of  the  different  mesh  sizes  used  in  the 
models,  coarse  grid  flows  must  be  proportioned  and  distributed  to  the  ap- 
propriate fine  grid  cells  in  order  for  flow  continuity  to  be  maintained. 

Flow  distribution  from  the  coarse  grid  to  the  fine  grid  system  is  made  on 
the  basis  of  the  relative  depths  of  the  fine  grid  cells  where  a particular 
coarse  grid  flow  is  to  be  specified.  For  example  if  the  flow  Qc  from  one 
coarse  grid  cell  is  to  be  input  to  the  fine  grid  sub-model  through  four 
cells  having  water  depths  of  di  = 1 , d2  =1 , d3  = 2 and  d4  = 2,  then  the 
corresponding  individual  fine  grid  flows  would  be  determined  as  Qi  = Q2  = 

(1/6)  Qc  and  Q3  = Q4  = (2/6)  Qc*  Both  the  temporal  interpolation  and  the 
spatial  distribution  of  coarse  grid  flows  are  carried  out  during  operation 
of  the  coarse  grid  model  with  the  input  flows  to  the  fine  grid  sub-model 
stored  on  magnetic  tape  at  the  end  of  each  run.  Table  8 summarizes  the 
linkage  between  the  two  sets  of  flows  in  the  models. 

Again  to  reflect  the  special  boundary  conditions  identified  in  Figure  24 
in  the  model,  the  computational  cell  flagging  scheme  shown  in  Figure  3 was 
used.  The  assigned  flags  used  in  the  fine  grid  sub-model  for  the  verifi- 
cation conditions  are  presented  in  Figure  25.  The  two-digit  cell  flags 
which  specify  in  the  model  appropriate  formulations  of  the  convective  acceler- 
ation cross-product  terms  are  shown  in  Figure  26.  The  values  of  Manning's  "n" 
bottom  roughness  coefficients  used  in  the  fine  grid  sub-model  were  determined 
in  a manner  consistent  with  that  described  earlier  for  the  coarse  grid  model. 
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Table  8.  Linkage  Between  Coarse  Grid  and  Fine  Grid  Boundary  Flows 


FINE  GRID  INPUT  FLOWS  COARSE  GRID  SOURCE  FLOWS 


Flow 

Number 

Cell 

I 

Indices 

J 

X/Y 

Component 

Cell 

I 

Indices 

J 

X/Y 

Component 

1 

9 

4 

Y 

5 

5 

Y 

2 

10 

4 

Y 

5 

5 

Y 

3 

11 

4 

Y 

5 

5 

Y 

4 

12 

4 

Y 

5 

5 

Y 

5 

13 

4 

Y 

6 

5 

Y 

6 

14 

4 

Y 

6 

5 

Y 

7 

15 

4 

Y 

6 

5 

Y 

8 

16 

4 

Y 

6 

5 

Y 

9 

17 

4 

Y 

7 

5 

Y 

10 

18 

4 

Y 

7 

5 

Y 

11 

19 

4 

Y 

7 

5 

Y 

12 

8 

5 

X 

4 

6 

X 

13 

8 

6 

X 

4 

6 

X 

14 

8 

7 

X 

4 

6 

X 

15 

8 

8 

X 

4 

6 

X 

16 

8 

9 

X 

4 

7 

X 

17 

8 

10 

X 

4 

7 

X 

18 

8 

11 

X 

4 

7 

X 

19 

8 

12 

X 

4 

7 

X 

20 

21 

13 

X 

7 

8 

X 

21 

21 

14 

X 

7 

8 

X 

22 

21 

15 

X 

7 

8 

X 

23 

1 

16 

X 

2 

9 

X 

24 

1 

17 

X 

2 

9 

X 

25 

1 

18 

X 

2 

9 

X 

26 

1 

19 

X 

2 

9 

X 

27 

32 

22 

X 

10 

10 

X 

28 

4 

23 

X 

3 

10 

X 

29 

32 

23 

X 

10 

10 

X 

30 

4 

24 

X 

3 

10 

X 

31 

32 

24 

X 

10 

10 

X 

32 

4 

25 

X 

3 

10 

X 

33 

20 

25 

X 

7 

11 

X 

34 

29 

25 

Y 

10 

10 

Y 

35 

31 

25 

Y 

10 

10 

Y 

36 

32 

25 

X 

10 

10 

X 

37 

32 

25 

Y 

10 

10 

Y 

38 

20 

26 

X 

7 

11 

X 

39 

20 

27 

X 

7 

11 

X 

40 

12 

28 

Y 

6 

11 

Y 

41 

13 

28 

Y 

6 

11 

Y 

42 

14 

28 

Y 

6 

11 

Y 

43 

15 

28 

Y 

6 

11 

Y 

44 

16 

28 

Y 

6 

11 

Y 

45 

20 

28 

X 

7 

11 

X 

46 

20 

29 

X 

7 

11 

X 

47 

20 

30 

X 

7 

11 

X 

Calibration  of  the  fine  grid  sub-model  followed  generally  the  same  pro- 
cedure as  that  used  for  the  coarse  grid  model.  Using  the  calibrated  coarse 
grid  model,  input  flows  to  the  fine  grid  sub-model  were  generated  for  each 
of  the  47  external  inflow  cells  identified  in  Figure  24  at  five  second  inter- 
vals over  the  entire  verification  period,  0400  to  2100  on  12  September  1969. 

For  the  first  run,  all  interior  flows  at  the  individual  cell  boundaries  were 
set  equal  to  zero,  but  initial  tidal  ele  ctions  were  interpolated  from  cor- 
responding coarse  grid  model  tides  to  assure  that  the  same  volume  of  water 
was  available  in  both  models  at  0400  hours.  For  subsequent  runs,  initial 
hydrodynamic  conditions  were  taken  as  those  computed  during  previous  operations 
of  the  fine  grid  sub-model  at  1650  hours  (one  complete  tidal  cycle).  In  this 
manner,  the  simulation  time  required  for  the  system  to  respond  to  a prescribed 
set  of  inputs  was  reduced. 

During  calibration  of  the  fine  grid  sub-model,  adjustments  were  made 
primarily  to  the  bottom  elevations  assigned  to  cells  in  the  major  channels 
and  to  the  barrier  elevations  and  discharge  coefficients.  Actually  very 
little  tuning  of  the  fine  grid  sub-model  was  required  to  adequately  simulate 
the  measured  velocities.  Final  results  from  the  model  are  shown  in  Figures 
27  to  41  in  terms  of  comparisons  between  computed  and  measured  velocities 
at  the  fifteen  stations  identified  in  Figure  4.  Velocities  in  the  model 
were  determined  by  dividing  the  flow  per  foot  of  width  across  selected  cell 
boundaries  by  the  average  water  depth  computed  from  water  levels  on  each 
side  of  the  boundaries.  This  procedure  was  necessary  because  of  the  manner 
in  which  variables  are  defined  on  the  discrete  elements  as  shown  in  Figure  2. 
Table  9 Is  a compilation  of  the  selected  cell  boundary  flows  which  were  used 
to  determine  velocities  in  the  model  for  comparison  with  prototype  data. 
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Table  9.  Cells  Used  in  the  Fine  Grid 
Sub-Model  for  Velocity  Verification 


Velocity 

Station 

Cell 

I 

Indices 

J 

Flow 

Component 

IN 

20 

10 

Y 

1C 

19 

10 

Y 

IS 

19 

10 

Y 

2N 

18 

17 

Y 

2C 

17 

17 

Y 

2S 

17 

17 

Y 

3N 

17 

24 

Y 

3C 

17 

24 

Y 

3S 

16 

24 

Y 

4E 

29 

22 

X 

4C 

29 

23 

X 

4M 

29 

24 

X 

5E 

5 

18 

X 

5C 

5 

18 

X 

5W 

5 

19 

X 

Although  the  simulation  of  prototype  tidal  conditions  with  the  fine  grid 
sub-model  was  considered  adequate  and  was  generally  consistent  with  that  of 
the  coarse  grid  model,  slight  differences  in  amplitude  and  phasing  were  noted 
at  some  gage  locations.  For  example  the  simulated  tide  at  Gage  5 was  very 
similar  In  both  shape  and  amplitude  to  the  prototype  tide,  but  it  occurred 
almost  an  hour  later.  At  Gage  1,  phasing  was  properly  simulated,  but  the 
computed  peak  tides  differed  from  prototype  measurements  by  as  much  as  0.4  of 
a foot.  Although  these  differences  are  not  considered  to  be  excessive  with 
regard  to  the  expected  accuracy  of  simulated  tides  and  particularly  because 
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of  the  good  agreement  between  computed  and  measured  velocities,  it  is  never- 
theless important  to  recognize  their  probable  causes  so  that  model  results 
can  be  interpreted  accordingly.  The  primary  reason  for  these  discrepancies 
lies  in  the  fact  that  the  coarse  grid  resolution  was  not  detailed  enough  to 
adequately  describe  all  of  the  inundated  areas  as  they  exist  in  the  prototype. 

At  some  locations  such  as  in  the  vicinity  of  Gage  5,  the  inundated  areas  could 
not  be  accounted  for  at  all  with  the  1200’  mesh  size,  even  though  the  ground 
elevation  at  other  inundated  cells  in  the  grid  were  adjusted  slightly  by 
appropriate  amounts  to  reflect  this  inundation  and  to  maintain  the  proper 
relation  between  depth  and  volume  in  the  model.  Even  with  these  adjusted 
elevations,  the  depth-volume  relations  in  the  two  models  could  not  be  made 
exactly  the  same  because  of  the  coarse  grid  resolution.  In  light  of  these 
considerations,  plus  the  fact  that  all  of  the  tide  gage  locations  are  not 
included  within  the  boundaries  of  the  fine  grid  sub-model,  it  is  felt  that 
tidal  simulations  with  the  coarse  grid  model  better  reflect  prototype  conditions. 

An  idea  of  the  general  velocity  pattern  in  the  immediate  vicinity  of  the 
inlet  can  be  determined  from  the  array  of  velocity  vectors  plotted  in  Figures 
42,  43,  44  and  45.  These  velocity  patterns  correspond  to  the  maximum  flood 
condition,  the  high  water  slack  condition,  the  maximum  ebb  condition,  and  the 
low  water  slack  condition,  respectively,  and  they  cover  an  area  bounded  by 
rows  5 and  24  and  columns  12  and  21  in  the  fine  grid.  Although  no  prototype 
measurements  are  available  with  which  to  verify  these  velocity  patterns,  they 
do  agree  with  anticipated  results,  particularly  with  regard  to  the  flow  direc- 
tions. 
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FLOOD  FLOW  CONDITIONS,  12  SEPTEMBER  1969 


SLACK  TIDE  CONDITION,  12  SEPTEMBER  1969 
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For  the  verification  period,  the  tidal  exchange  in  the  prototype  was 
determined  to  be  641  million  cubic  feet  from  the  velocity  measurements  made 
at  Range  2 and  the  tidal  history  recorded  at  Gage  1.  From  the  fine  grid  sub- 
model results  at  this  same  velocity  range,  an  exchange  of  678  million  cubic 
feet  was  computed.  Both  of  these  exchanges  were  determined  for  the  flood 
flow  period  1400  to  2100  hours  on  12  September  1969.  The  difference  between 
the  two  values  is  approximately  five  percent,  and  an  error  of  this  magnitude 
is  not  considered  excessive,  especially  in  light  of  the  approximations  made 
in  the  model  and  the  accuracy  with  which  field  measurements  can  be  made. 

The  degree  to  which  the  final  results  from  both  the  fine  and  coarse 
grid  models  compare  with  actual  prototype  measurements  directly  reflects 
the  level  of  effort  available  in  this  study  for  calibration  of  the  models. 

It  is  possible  that  still  further  adjustment  and  tuning  of  the  models  could 
provide  better  reproduction  of  the  prototype  tides  and  velocities,  however, 
it  Is  believed  that  these  results  adequately  demonstrate  the  simulation 
capabilities  of  the  models  and  provide  some  measure  of  the  level  of  confi- 
dence with  which  the  computed  results  can  be  used.  Certainly  finer  resolution 
in  both  the  models  would  allow  better  descriptions  of  physical  detail,  and 
undoubtably  better  simulation  of  prototype  conditions.  However,  it  is  felt 
that  the  mesh  sizes  selected  are  the  absolute  minimums  which  could  have  been 
used  in  this  study  without  model  operation  costs  becoming  excessive. 
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IV.  PRE-  AND  POST-PROJECT  SIMULATIONS 


Model  Modifications 

With  the  calibration  and  verification  of  the  models  for  the  12  Sep- 
tember 1969  conditions  satisfactorily  completed,  operation  of  the  models 
was  undertaken  to  simulate  the  inlet  hydrodynamics  for  the  following  two 
conditions: 


1.  Pre-project  conditions  (without  jetty ) aa  of  November  1964,  and 

2.  Post-project  conditions  (with  jetty)  as  of  June  1967. 

Each  of  these  inlet  configurations  was  investigated  using  both  mean  and 
spring  tides  to  excite  the  models.  Figure  46  illustrates  these  tides  for 
one  semi-diurnal  tidal  cycle. 

Physiographic  information  for  each  condition  was  provided  by  means  of 
1"  = 500'  hydrographic  surveys  covering  the  imnedlate  inlet  area.  The  spatial 
limit  of  these  surveys  is  shown  in  Figure  47  relative  to  the  model  grids. 

The  Information  provided  did  not  extend  over  either  of  the  entire  grid  areas, 
and  therefore,  for  simulation  of  both  June  1967  and  November  1964  conditions, 
the  bottom  and  ground  elevations  and  water-land  boundaries  outside  of  the 
survey  limits  were  not  changed  in  the  models  from  those  used  In  the  verifi- 
cation process.  Alterations  were  made  only  to  Interior  cells. 

Five-foot  bottom  contours  were  provided  for  each  of  the  operational 
conditions  on  the  hydrographic  surveys  together  with  MLW  and  MHW  lines. 
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While  no  ground  elevations  were  specified,  it  is  unlikely  that  the  surface 
topography  of  1964  and  1967  was  significantly  different  from  that  described 
on  the  1969  hydrographic  chart.  For  the  June  1967  condition,  the  location 
of  the  jetty  and  the  weir  section  were  also  provided.  The  specific  changes 
that  were  made  to  each  of  the  model  grids  and  the  attendent  input  data  alter- 
ations for  simulation  of  the  two  operational  conditions  are  described  below. 

November  1964  Condition.  Of  the  two  operational  conditions,  the 
November  1964  inlet  geometry  required  the  most  modifications  to  the  grid 
networks  used  in  verifying  the  models.  The  major  change  was  the  removal  of 
the  jetty  and  the  inclusion  In  the  fine  grid  of  several  additional 
water  cells  in  the  immediate  jetty  area.  New  bottom  and  ground  elevations 
were  assigned  to  cells  in  both  grids  as  required.  Also  new  computational 
cell  flags  and  convective  acceleration  cell  flags  were  specified  where  changes 
in  cell  boundary  conditions  were  made.  Without  the  jetty  in  place,  the  normal 
flow  of  water  was  allowed  across  the  original  jetty  alignment,  and  cell  flags 
were  modified  to  reflect  these  conditions.  In  the  fine  grid  sub-model,  ten 
additional  external  inflow  boundaries  were  created  to  account  for  these  flows 
into  the  system. 

Where  changes  in  bottom  depths  resulted  in  the  creation  or  removal  of 
inundated  cells,  flags  were  changed  to  handle  these  new  flow  conditions. 
Submerged  barrier  elevations  were  made  consistent  with  surrounding  bottom 
topography. 

Because  of  the  limited  hydrographic  survey  information  Included  in 
Figure  47,  the  changes  in  the  coarse  grid  configuration,  Figure  8,  affected 
only  cells  between  rows  5 and  11  and  columns  4 and  8.  Similarly,  the  fine 
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grid  changes  were  made  only  to  cells  between  rows  4 and  24  and  columns  9 and 
22,  Figure  22.  The  new  arrays  of  bottom  and  ground  elevations,  computational 
cell  flags,  and  convective  cell  flags  for  both  models  are  presented  in  Figures 
48  through  53  for  the  November  1964  conditions.  Table  10  describes  the  linkage 
between  the  two  models  for  the  ten  additional  external  inflow  boundary  cells 
created  in  the  fine  grid  sub-model. 

Table  10.  Linkage  Between  Additional  Coarse  Grid  and  Fine  Grid  Boundary 
Flows  for  November  1964  Conditions 


FINE  GRID  INPUT  FLOWS  COARSE  GRID  SOURCE  FLOWS 


Flow 

Number 

Cell 

I 

Indices 

J 

X/Y 

Component 

Cell 

I 

Indices 

J 

X/Y 

Component 

12 

20 

4 

Y 

7 

5 

Y 

13 

21 

4 

Y 

7 

5 

Y 

15 

21 

5 

X 

7 

6 

X 

17 

21 

6 

X 

7 

6 

X 

19 

21 

7 

X 

7 

6 

X 

21 

21 

8 

X 

7 

6 

X 

23 

21 

9 

X 

7 

7 

X 

25 

21 

10 

X 

7 

7 

X 

27 

21 

11 

X 

7 

7 

X 

29 

21 

12 

X 

7 

7 

X 

June  1967  Condition.  Because  the  inlet  geometry  during  this  period  was 
quite  similar  to  that  used  for  calibration  of  the  models,  only  a few  modifications 
were  required  to  describe  the  physiography  of  this  condition.  Of  the  changes 
made  all  were  required  because  of  the  different  bottom  topography,  and  conse- 
quently external  boundary  conditions  in  both  the  coarse  and  the  fine  grid  models 
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Figure  48 

MSL  BOTTOM  AND  GROUND 
ELEVATIONS  USED  IN  COARSE 
GRID  MODEL  FOR  NOVEMBER 
1964  CONDITIONS 

4 5 6 7 8 

X 

Figure  49 

COMPUTATIONAL  CELL  FLAGS  USED 
IN  COARSE  GRID  MODEL  FOR 
NOVEMBER  1964  CONDITIONS 


4 5 6 7 8 

X 


Figure  50 

CONVECTIVE  ACCELERATION  CELL 
FLAGS  USED  IN  COARSE  GRID 
MODEL  FOR  NOVEMBER  1964 
CONDITIONS 

4 5 6 7 8 
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remained  the  same.  Because  bottom  elevation  changes  caused  some  cells  in  the 
models  to  become  periodically  inundated,  computational  cell  flags  had  to  be 
modified  accordingly.  Also,  a few  barrier  elevations  were  changed  to  reflect 
corresponding  changes  in  bottom  elevations. 

Using  the  same  localized  grid  areas  as  those  described  earlier  for  the 
November  1964  model  modifications,  the  new  bottom  and  ground  elevations, 
computational  cell  flags,  and  convective  cell  flags  required  to  simulate 
the  June  1967  conditions  are  presented  in  Figures  54  through  59  for  both 
the  coarse  and  fine  grid  models. 

Tidal  Excitation  in  Coarse  Grid  Model 


No  modifications  to  the  tidal  excitation  boundaries  in  the  coarse  grid 
model  were  required  by  either  of  the  two  operational  conditions.  Both  utilized 
the  same  tidal  input  cells  designated  in  Figure  4;  however,  new  tidal  input 
curves  had  to  be  constructed  from  the  basic  mean  and  spring  ocean  tides  given 
in  Figure  46.  To  determine  the  three  inland  exciting  tides  corresponding  to 
the  mean  and  spring  tides,  the  tidal  constants  developed  during  calibration 
of  the  model  and  presented  in  Table  4 were  used  in  eq.  (32)  with  the  Fourier 
coefficients  and  phase  angles  shown  in  Table  11.  The  resulting  exciting 
tides  used  for  simulation  of  both  the  November  1964  and  the  June  1967  con- 
ditions are  presented  in  Figures  60  and  61  for  the  mean  and  spring  tides, 
respectively.  These  tides  were  assigned  to  appropriate  boundary  cells  according 
to  the  tidal  identification  numbers  given  in  Table  3. 
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MEAN  EXCITING  TIDES  USED  IN  COARSE  GRID  MODEL 


Table  11.  B1 -Modal  Fourier  Constants 
Used  to  Approximate  the  Mean  and  Spring  Tides 


Parameter 

Mean  Tide 

Spring  Tide 

Ao 

1 .88200 

1 .85040 

Co  i 

1.93147 

2.28753 

$0  1 

-1.54206 

-1.54514 

Co2 

0.03943 

0.03028 

$02 

-1.17606 

-1.38140 

Operational  Results 

Once  the  necessary  changes  had  been  made  in  the  models  to  accommodate  the 
physiographic  conditions  of  the  1964  and  1967  operational  periods,  simulation 
of  the  hydrodynamics  was  carried  out  using  the  same  procedure  as  that  employed 
in  the  verification  simulations.  For  each  inlet  geometry  and  prescribed  ex- 
citing tide,  the  coarse  grid  model  was  operated  to  obtainthe  necessary  boundary 
input  flows  for  the  fine  grid  sub-model.  With  the  initial  hydrodynamics  set 
equal  to  zero,  two  complete  tidal  cycles  (=25  hours)  were  simulated  with  the 
flows  and  tides  from  the  second  cycle  used  to  determine  the  Inputs  to  the  fine 
grid  sub-model.  Again  the  coarse  grid  flows  were  Interpolated  temporally  and 
distributed  spatially  to  provide  the  proper  linkage  with  the  fine  grid  flows. 
Time  steps  in  both  models  were  the  same  as  those  used  in  verification.  With 
the  fine  grid  exciting  flows  stored  on  magnetic  tape,  the  fine  grid  sub-model 
was  operated  for  two  tidal  cycles  beginning  with  initial  hydrodynamics  set 
equal  to  zero.  Results  from  the  second  cycle  were  used  for  presentation. 
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For  both  of  the  operational  periods,  the  average  wind  condition  for  each 
of  the  months  was  included  in  the  models  for  all  of  the  simulations.  The  wind 
data  was  obtained  from  the  NOAA  Environmental  Data  Service  Division  in  Asheville, 
North  Carolina,  and  the  wind  conditions  measured  at  the  Wilmington,  North 
Carolina,  first  order  station  were  used  in  the  models.  The  wind  speed  and 
direction  were  held  constant  throughout  the  entire  simulation  period. 

The  results  from  the  simulations  for  each  of  the  two  inlet  conditions 
with  the  mean  and  spring  exciting  tides  are  presented  in  the  Appendix.  Tidal 
elevations  computed  with  the  models  are  shown  for  each  of  the  six  tide 
gage  locations  in  the  inlet  area.  The  same  coarse  grid  cells  identified  in 
Table  6 were  used  to  obtain  the  computed  tides  at  the  gage  locations.  Gage  0 
results  are  not  presented  since  it  has  been  assumed  that  the  prescribed  off- 
shore exciting  tides  (Figure  46)  correspond  to  those  which  would  occur  at  the 
Gage  0 location.  The  computed  velocities  are  also  shown  at  the  same  stations 
identified  in  Figure  4 which  were  used  for  model  verification.  Because  the 
three  stations  at  several  of  the  velocity  ranges  were  too  closely  spaced  to 
allow  each  to  be  described  with  a separate  cell  in  the  fine  grid,  the  velocity 
simulated  for  a single  grid  cell  in  some  instances  had  to  be  used  as  the 
velocity  which  occurs  at  two  adjacent  stations.  For  this  reason,  some  of  the 
velocity  plots  in  the  Appendix  correspond  to  two  velocity  stations,  and  they 
are  labeled  accordingly. 

Using  the  mean  and  spring  tidal  flows  computed  across  Range  2 in  the 
fine  grid  sub-model,  CELLS  (17,17),  (18,17),  and  (19,17),  the  total  exchange 
during  the  flood  period  was  determined  for  both  the  November  1964  and  June 
1967  inlet  geometries.  For  the  1964  condition,  the  mean  tidal  exchange  was 
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determined  to  be  489  million  cubic  feet  and  the  spring  tidal  exchange  was 
computed  as  565  million  cubic  feet.  Similarly  for  the  1967  inlet  geometry 
with  the  jetty  in  place,  mean  and  spring  tidal  exchanges  were  calculated 
to  be  400  and  451  million  cubic  feet,  respectively.  These  smaller  values 
are  consistent  with  those  anticipated,  and  they  reflect  the  changes  made  to 
ground  elevations  in  each  of  the  models  for  the  two  operational  conditions, 
as  well  as  the  reduced  influence  of  ocean  exciting  tides  with  the  jetty 
in  place. 

Generalized  velocity  patterns  at  the  maximum  flood  flow  condition  for 
both  the  mean  and  spring  tides  are  shown  in  Figures  62  through  65  for  each 
of  the  inlet  geometries.  Only  slight  differences  are  discernible  between 
the  mean  and  spring  tidal  conditions  for  the  same  Inlet  geometry,  however 
as  would  be  expected,  there  are  very  pronounced  changes  In  the  velocities 
for  the  simulations  with  and  without  the  jetty  In  place.  For  the  pre-project 
November  1964  condition,  long-shore  velocities  are,  of  course,  somewhat 
higher  than  those  computed  with  the  jetty  and  the  result  of  these  higher 
velocities  is  directly  reflected  in  the  exchange  values  presented  above  for 
the  two  inlet  configurations. 
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FLOW  CONDITION  FOR  SPRING  TIDE,  NOVEMBER  1964 
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FLOW  CONDITION  ^OR  SPRING  TIDE,  JUNE  1967 


V.  SUMMARY  AND  RECOMMENDATIONS 


A spatially  two-dimensional  tidal  hydrodynamics  model,  HYDTID,  has 
been  applied  to  the  Masonboro  Inlet  tidal  system  and  calibrated  and  verified 
for  prototype  conditions  of  12  September  1969.  Model  verification  for  this 
period  consisted  of  reproducing  semi-diurnal  tides  at  six  gage  locations, 
velocities  at  three  stations  along  each  of  five  ranges,  and  tidal  exchange 
at  the  mouth  of  the  inlet. 

In  structuring  the  model,  a dual-model  approach  has  been  employed 
whereby  a coarse  grid  model  covering  the  entire  inlet  area  is  used  to  gener- 
ate the  exciting  flows  for  a fine  grid  sub-model  of  the  immediate  inlet  area. 
The  course  grid  model  utilizes  a 1200  foot  mesh  size  to  describe  the  inlet 
geometry  and  physical  characteristics  of  the  system  whereas  the  fine  grid 
sub-model  is  based  on  a 300  foot  grid  spacing.  Linkage  between  the  two 
models  is  provided  along  common  cell  boundaries  where  the  coarse  grid  flows 
are  temporally  interpolated  and  spatially  distributed  to  obtain  the  input 
flows  to  the  fine  grid  sub-model.  Excitation  of  the  coarse  grid  model  is 
accomplished  by  prescribing  predetermined  tidal  elevations  around  the  periph- 
ery of  the  model . 

Calibration  and  verification  was  first  undertaken  for  the  coarse  grid 
model  with  primary  emphasis  placed  on  tidal  simulation  and  the  general  distri- 
bution of  flows  within  the  major  channels  and  passes  in  the  inlet  area.  The 
fine  grid  sub-model  was  then  calibrated  to  simulate  the  localized  velocities 
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at  each  of  the  stations  along  the  velocity  ranges.  An  important  consideration 
in  calibrating  each  of  the  models  was  that  the  relation  between  water  depth 
and  volume  be  the  same  in  both  of  the  models  for  that  portion  of  the  prototype 
included  in  the  fine  grid  sub-model. 

With  both  of  the  models  satisfactorily  calibrated  and  verified,  operational 
runs  were  carried  out  to  simulate  the  inlet  hydrodynamics  under  mean  and  spring 
tidal  conditions  for  two  different  inlet  geometries.  One  series  of  runs  corre- 
sponded to  the  November  1964  pre-project  condition  without  the  jetty,  and  the 
other  was  for  the  modified  inlet  configuration  of  June  1967  with  the  jetty  in 
place.  Alterations  were  made  to  each  of  the  models  to  reflect  the  new  inlet 
conditions  and  simulations  were  made  for  two  complete  tidal  cycles.  Results 
were  compared  at  the  same  tide  gage  locations  and  velocity  stations  used  for 
verifying  the  models. 

It  is  believed  that  the  conformity  of  the  simulations  obtained  in  this 
study  are  consistent  with  the  numerical  techniques  employed  in  the  model  and 
the  accuracy  with  which  field  measurements  can  be  made.  While  it  is  possible 
to  obtain  even  better  agreement  between  the  computed  and  measured  values  by 
making  further  model  adjustments,  small  incremental  improvements  can  be 
achieved  only  at  the  expense  of  additional  calibration  runs.  Such  extended 
verification  may  not  be  justified  without  at  least  some  supplemental  flow 
information  at  some  other  locations  near  the  inlet.  Also  it  would  be  desirable 
to  extend  the  model  verification  to  other  tidal  conditions.  With  respect  to 
operation  of  the  Masonboro  Inlet  models,  the  following  specific  recomnendations 
should  be  considered. 
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1.  As  additional  prototype  information  becomes  available  the  models 
should  be  applied  arid  adjusted  as  required  to  extend  their  range 
of  verification. 

2.  The  correlations  between  the  offshore  tide  and  those  prescribed 
at  the  inland  artificial  tidal-storage  boundaries  should  be 
tested  for  other  tidal  conditions. 

3.  Although  the  coarse  and  fine  grid  models  in  their  present  form 
use  time  steps  of  twenty  and  five  seconds , respectively , it  is 
suggested  that  increments  of  twenty-four  and  six  seconds  be 
used  if  it  is  desired  to  reduce  computational  run  time. 

For  future  applications  of  HYDTID  to  other  tidal  inlets,  consideration 
should  be  given  to  the  following. 

1.  In  order  to  facilitate  the  application  of  HYDTID,  concerted 
effort  should  be  made  to  collect  the  information  necessary  to 
establish  model  boundary  conditions  rather  to  develop  them  in 
the  calibration  process. 

2.  When  employing  the  dual-model  approach  it  would  be  desirable 
to  incorporate  in  the  program  an  internal  check  on  the  depth- 
volume  relations  of  the  two  models. 

3.  The  HYDTID  program  could  be  refined  by: 

a.  Elimination  of  those  subroutines  not  required  for  a 
specific  simulation  of  inlet  hydrodynamics,  and 

b.  Further  generalization  of  the  program  to  provide  for 
direct  input  of  all  required  model  data. 

The  Masonboro  Inlet  problem  is  certainly  one  of  the  more  complex  appli- 
cations of  HYDTID.  The  major  difficulty  encountered  was  providing  the  necessary 
spatial  resolution  to  describe  the  local  detail  in  a cost-effective  manner. 

While  the  methodology  employed  proved  workable,  it  is  probable  that  a more 
efficient  procedure  could  be  developed  within  the  HYDTID  framework.  Such 
procedure  might  involve  a variable  grid  or  a series  of  non-overlapping  sub- 
models with  different  mesh  sizes  systematically  linked  through  the  hydrodynamics. 
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Appendix  A 

SIMULATED  RESULTS  FOR  NOVEMBER  1964  INLET  GEOMETRY 
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SIMULATED  RESULTS  FOR  JUNE  1967  INLET  GEOMETRY 
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